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ABSTRACT 

Using the extensive multi-wavelength data in the GOODS-North field, including our ground-based 
rest-frame UV spectroscopy and near-IR imaging, we construct and draw comparisons between samples 
of optical and near-IR selected star-forming and passively evolving galaxies at redshifts 1.4 ^ z < 2.6. 
We find overlap at the 70 — 80% level in samples of 2; '--^ 2 star-forming galaxies selected by their 
optical (UnGTZ) and near-IR (BzK) colors when subjected to common K-hand limits. Deep Chandra 
data indicate a ~ 25% AGN fraction among near-IR selected objects, much of which occurs among 
near-IR bright objects {K^ < 20; Vega). Using X-rays as a proxy for bolometric star formation 
rate (SFR) and stacking the X-ray emission for the remaining (non-AGN) galaxies, we find the SFR 
distributions of UnGTZ, BzK, and J ~ > 2.3 galaxies (i.e., Distant Red Galaxies; DRGs) are 
very similar as a function of Kg, with A's < 20 galaxies having (SFR) ~ 120 Mq yr~^, a factor of 
two to three higher than those with > 20.5. The absence of X-ray emission from the reddest 
DRGs and BzK galaxies with (z — A')ab ^ 3 indicates they must have declining star formation 
histories to explain their red colors and low SFRs. While the M/L ratio of passively-evolving galaxies 
may be larger on average, the Spitzer/IRAC data indicate that their inferred stellar masses do not 
exceed the range spanned by optically selected galaxies, suggesting that the disparity in current 
SFR may not indicate a fundamental difference between optical and near-IR selected massive galaxies 
(M* > 10^^ M0). We consider the contribution of optical, near-IR, and submillimeter-selected galaxies 
to the star formation rate density (SFRD) at z '--^ 2, taking into account sample overlap. The SFRD 
in the interval 1.4 < z < 2.6 of UnGTZ and BzK galaxies to K^ = 22, and DRGs to K^ = 21 is 
- 0.10 ± 0.02 Mq yr~i Mpc"^. Optically-selected galaxies to 7^ = 25.5 and K^ = 22.0 account for 
~ 70% of this total. Greater than 80% of radio-selected submillimeter galaxies to S'sso/im ~ 4 mJy 
with redshifts 1.4 < z < 2.6 satisfy either one or more of the BX/BM, BzK, and DRG criteria. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high redshift — galaxies: 
starburst — infrared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

A number of surveys have been developed to select 
galaxies at z ~ 2, determine their bolometric star for- 
mation rates (SFRs), and compare with other multi- 
wavelength studies to form a census of the total star for- 
mation rate density (SFRD) at z ~ 2 fe.g.. 'S teidel et aO 
|1]04; Rubin ct al. 2004; Daddi et al. 2004b). A paral- 
lel line of study has been to compare optical and near- 
IR selected galaxies that are the plausible progenitors 
of the local population of passively evolving massive 
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ifornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
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galaxies. However, biases inherent in surveys that se- 
lect galaxies based on their star for matio n act i vity (e.g.. 
'Steidel et al."2004^ and stellar mass fe.g.. iCimatti et all 
2002b; Glazebrook et al. 2004) can complicate such com- 
parisons. Only with an accurate knowledge of the over- 
lap between these samples can we begin to address the 
associations between galaxies selected in different ways, 
their mutual contribution to the SFRD at z ~ 2, and the 
prevalence and properties of passively evolving and mas- 
sive galaxies at high redshift. Quantifying this overlap 
between optical and near-IR surveys is a primary goal of 
this paper. 

In practice, optical surveys are designed to efficiently 
select galaxies with a specific range of properties. The 
imaging required for optical selection is generally a small 
fraction of the time required for near-IR imaging, and can 
cover much larger areas within that time. In contrast, 
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near-IR surveys sample galaxies over a wider baseline in 
wavelength than optical surveys, and can include galax- 
ies relevant to studying both the star formation rate and 
stellar mass densities at high redshift. However, in order 
to achieve a depth similar (and area comparable) to that 
of optical surveys, near-IR selection requires extremely 
deep imaging and can be quite expensive in terms of tele- 
scope time due to the relatively small size of IR arrays 
compared to CCDs. Furthermore, the "color" of the ter- 
restrial background for imaging is {B—Ks)ab — 7 magni- 
tudes, much redder than all but the most extreme z ^ 2 
galaxies. Once selected, of course, such extreme galax- 
ies then require heroic efforts to obtain spectra, whereas 
optical selection, particularly at redshifts where key fea- 
tures fall shortward of the bright OH emission "forest" , 
virtually guarantees that one can obtain a spectroscopic 
redshift with a modest investment of 8-lOm telescope 
time and a spectrograph with reasonably high through- 
put. As we show below, optical and near-IR surveys 
complement each other in a way that is necessary for ob- 
taining a reasonably complete census of galaxies at high 
redshift. 

The SFRs of z ~ 2 galaxies are typically estimated 
by employing locally-calibrated relations between emis- 
sion at which the galaxies can be easily detected (e.g., 
UV, Ha) and their FIR emission. The X-ray luminos- 
ity of local non-active galaxies results primarily from 
high mas s X-ray binaries, supernovae, and diffuse hot 
gas re.g.. lGrimm et al.ll2(m iStrickland et al.ll2f)0l : all 
of these sources of X-ray emission are related to the star 
formation activity on timescales of < 100 Myr. Observa- 
tions of galaxies in the local universe show a tight corre- 
lation between X-ray and FIR luminosity, prompting the 
use o f X-ray emission as an SFR indicator l)R,analli et al.l 
l2003fl . This correlation between X-ray emission and SFR 
applies to galaxies with a very large range in SFRs, from 
^ 0.1—1000 Mq yr~^. Stacking analyses at X-ray and ra- 
dio wavelengths, and comparison with UV emission, indi- 
cate that the local SFR relations appear to give compara- 
ble estimates of the instantaneous SFRs of galaxies after 
assuming continuous star formation models and correct- 
ing for dust Ce.g.. iReddv fc Steidell I2004t iNandra et all 
l2nr)l ISeibert et a,lJl2002D ! 

Two surveys designed to select massive galaxies at red- 
shifts 1.4 S z < 2.5 and passively-evolving (PE) galaxies 
at redshifts z > 2, respectively, are the K20 and FIRES 
surveys. The K20 and FIRES selection criteria were de- 
veloped to take advantage of the sensitivity of rest-frame 
optical light and color to stellar mass and the strength of 
the Balmer break, respecti vely, for z ~ 2 galaxies (e.g., 
IHimatti ct al. 2002a; Fran x et al J 1200^1) . The Gemini 
Deep Deep Survey (GDDS) extends this near-IR tech- 
nique to target massive galaxies at sli ghtly lower redshifts 
(0-8 <z< 2.0: IAbraham et al.ll2004tl . 

X-ray stacking analyses of the brightest galaxies in the 
K20 and FIRES surveys indicate an average SFR a fac- 
tor of 4 to 5 times larger than fo r optically -selected z ~ 2 
galaxies ( Daddi et ah,.2004a: .R.ubin et al.li2004 'l , inviting 
the conclusion that optical selection misses a large frac- 
tion of the star formation density at high redshift. While 
it is certainly true that optical surveys miss some fraction 
of the SFRD, the past quoted difference in the average 
SFRs of galaxies selected optically and in the near-IR 
disappears once the galaxies are subjected to a common 



near-IR magnitude limit, as we show below. 

We have recently concluded a campaign to obtain deep 
near-IR im aging for fields in the z ~ 2 optical survey 
llSteidel et al. 2004), allowing for a direct comparison of 
optical and near-IR selected galaxies. One result of this 
comparison is that Kg < 20 (Vega) optically-selected 
galaxies show similar space densities, stellar masses, and 
metallicitics as ifg-bright galaxies in jiear-IR s amples 
Iphaolcv ct al. 2004). More recently, Adclbcrecr et alJ 
(2005) show that the correlation lengths for Xg-bright 
galaxies among optical and near-IR samples are similar, 
suggesting an overlap between the two sets of galaxies, 
both of which plausibly host the progenitors of massive 
elliptical galaxies in the local universe. These results sug- 
gest that near-IR bright galaxies have similar properties 
regardless of the method used to select them. 

In this paper, we extend these results by examining 
the color distributions and X-ray properties of near-IR 
and optically selected galaxies at z ~ 2 in the GOODS- 
North field ( Giava lisco et al.ll2004D . The field is well- 
suited for this analysis given the wealth of complemen- 
tary data available, including Chandra/ ji-ray, ground- 
based optical and near-IR, and Spitzer/IKAC imaging. 
Multi-wavelength data in a single field are particularly 
useful in that we can use a common method for extract- 
ing photometry that is not subject to the biases that may 
exist when comparing galaxies in different fields whose 
fluxes are derived in different ways. The addition of our 
rest-frame UV spectroscopic data in the GOODS-N field 
provides for a more detailed analysis than otherwise pos- 
sible of the properties of galaxies as a function of selec- 
tion technique. Furthermore, the GOODS-N field coin- 
cides with the Chandra Deep Field North (CDF-N) re- 
gion which have the de epest (2 Ms) X-ray data available 
|! Alex ander et al.ll2003(l . The X-ray data allow for an in- 
dependent estimate of bolometric SFRs and the available 
depth allows more leeway in stacking smaller numbers of 
sources to obtain a statistical detection, as well as identi- 
fying AGN to a lower luminosity threshold than possible 
in other fields that have shallower X-ray data. 

The outline of the paper is as follows. In §2, we de- 
scribe the optical, near-IR, X-ray, and IRAC data and 
present the optical and near-IR selection criteria and X- 
ray stacking method. Color distributions, direct X-ray 
detections, and stacked results are examined in §3. In §4, 
we discuss the SFR distributions of optical and near-IR 
selected z ~ 2 galaxies and their relative contributions to 
the SFRD, and the presence of a passively evolving pop- 
ulation of galaxies. A flat ACDM cosmology is assumed 
with Ho = 70 km s"^ Mpc"^ and f^A = 0.7. 

2. DATA AND SAMPLE SELECTION 

2.1. Imaging 

Optical UnGTZ images in the GOODS-North flcld were 
obtained in 2002 and 2003 April under photometric con- 
ditions using the KPNO and Keck I telescopes. The 
KPNO/MOSAIC [/-band image was obtained from the 
GOODS team (PI: Giavalisc o) and was transfor med to 
reflect C/„ magnitudes (e.s.. iSteidel et al.ll200l . The 
Keck I G and TZ band images were taken by us wi th the 
Low Resolution Imaging Spectrograph (LRIS; Oke et all 
ISteidel et al.' 2004), and were oriented to provide 
the maximum overlap with the GOODS ACS and Spitzer 
survey region. The images cover 11' x 15' with FWHM 
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Table 1. Interloper Contamination of the BX/BM Sample 



Ks Range 


A^phot" 


Af t> 
i 'spec 




f , d 
Jz<l 


< 20.0 


61 


18 


7 


0.39 


20.0 < Ks< 20.5 


58 


23 


2 


0.09 


20.5 < i^s < 21.0 


82 


30 


3 


0.10 


21.0 < i^s < 21.5 


101 


32 


2 


0.06 


21.5 < Ks < 22.0 


141 


29 


3 


0.10 



''Number of photometric BX/BM candidates. 
''Number of candidates with spectroscopic redshifts. 
'^Number of interlopers. 
"^Interloper fraction. 



~ O'.'T to a depth of 7?. ~ 27.5 (3 a). Image reduction and 
p hotometry wer e done following the procedures described 
in iSteidel et al .> ("200?). We obtained deep B-band im- 
ages of the GOODS-N fiel d from a public distribution of 
Subaru data (fCapak et al. 2004). The deep z-band data 
are acquired from the public distributio n of the HST Ad- 
vance d Camera for Surveys (AGS) data ijGiavalisco et al.l 
I2004D . The B and z band data have 5 a depth of 26.9 and 
27.4 mag measured in 3 " and 0'.'2 diameter apertures, 
respectively. The Ks and J imaging was accomplished 
with the Wide Field Infrared Gamera (WI RC) on the 
Palomar Hale 5 m telescope l|Wilson et al.l 12003'). pro- 
viding 8.'7 X 8.'7 coverage in the central portion of the 
GOODS-N field. The near-IR images cover - 43% of the 
optical image. The images had FWHM ^ I'.'O under pho- 
tometric conditions and 3 a sensitivity limits of ~ 22.6 
and 24.1 mag in the Kg and J bands, re spectively. 
The n ear-IR data are described in detail by lErb et aTl 
l|200l . The total area studied in the subsequent analy- 
sis is ~ 72.3 arcmin^. 

The procedure s for source d e tectio n and photometry 
are described in ISteidel et all l|2003t) . Briefly, UnGTl 
magnitudes were calculated assuming isophotal aper- 
tures that were adjusted to the 7?,-band flux profiles. 
Source detection was done at Ks-hand. BzK and J mag- 
nitudes are computed assuming the isophotal apertures 
adjusted to the i^g-band flux profiles, unless the 7?.-band 
isophotes gave a more significant Kg detection. In the 
analysis to follow, "i^s" and J magnitudes are in Vega 
units. We use the conversion Kab = + 1-82. All other 
magnitudes are in AB units. 

Fully reduced Spitzer/IKKC mosaics of the GOODS- 
North field were made public in the first data release of 
the GOODS Legacy project (PI: Dickinson). The IRAG 
data overlap completely with our i^Tg-band image, but 
currently only two channels (either 3.6 /.tm and 5.8 /xm, 
or 4.5 /im and 8.0 ^m) are available over most of the 
image. A small area of overlap has coverage in all four 
channels. The images are deep enough that source con- 
fusion is an issue. We have mitigated the effects of con- 
fusion noise by employing the higher spatial resolution 
Xg-band data to constrain source positions and deblend 
confused IRAG sources. We perf ormed PSF photomet ry 
using the procedure described in lShaDlev et al.l l|200,'if) . 

2.2. Selection Criteria 
Optical Selection of Star-Forming Galaxies 



We have optically-selected z ~ 2 galaxies in the 
GOODS-N field based on their observed UnGTZ colors 
ijAdelberger et al.ll2004t Isteidel et al.ll2004|) to a hmiting 
magnitude oiTZ = 25.5. The selection criteria aim to se- 
lect actively star forming galaxies at 2 ~ 2 with the same 
range in UV p roperties and extinction as LBGs at z ~ 3 
(jSteidel et alJ l2b03). "BX" galaxies are selected to be at 
redshifts 2.0 < z < 2.6 using the following criteria: 

G-7^>-0.2 
Un-G>G-n + 0.2 

G-7^<0.2(C/„-G)^- 0.4 
[/„-G<G-7^ + 1.0, (1) 

and "BM" objects are selected to be at redshifts 1.5 < 
z < 2.0 using the following criteria: 

G-7^>-0.2 
^7„-G>G-7^-0.1 

G-7^<0.2([/„-G)-l-0.4 
Un-G<G-n + 0.2 (2) 

l|Adelberger et al.l l2004: Stcidcl ct al."2004V For subse- 
quent analysis, we will refer to BX and BM objects as 
those which are optically-, or "BX/BM"-, selected. 

Optical color selection of z ^ 2 galaxies in the 11' by 
15' area of the GOODS-North field yielded 1360 BX and 
BM candidates, of which 620 lie in the region where we 
have complementary J and i^-band data fS I2.1|l . and 199 
have Kg < 21.0. FoUowup spectroscopy with the blue 
channel of the Low Resolution Imaging Spectrograph 
(LRIS-B) yielded 147 redshifts for objects with i^Tg-band 
data (248 redshifts over the entire optical field). Of these 
147 objects with redshifts and iiTg-band data, 129 have 
z > 1, and 60 have z > 1 and Ks < 21. The mean red- 
shift of the 60 BX/BM objects is (z) = 1.99 ± 0.36. The 
spectroscopic interloper fractions in the BX/BM sample 
are summarized in Table The BX and BM selection 
functions (shown as shaded distributions in Figure ^ 
have distributions with mean reds hifts (z) = 2.2 ± 0. 3 
and (z) = 1.7 ±0.3, respectively (|Steidel et al.ll2004() . 
and the combination of these two samples comprise our 
BX/BM-selected z ~ 2 sample. In the analysis to follow, 
we designate an interloper as any object with z < 1. 

Near-IR Selection of Star-Forming Galaxies 

The near-IR properties of galaxies can be used both to 
target star forming galaxies and to identify those with 
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BX/BM: (z)=1.99±0.36 ^ 
BzK/SF: (z) = 2.01±0.33 ^ 
BzK/PE: (z>~1.70±0.18 
DRG: (z) = 2.48±0.38 
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z 



3.5 



Fig. 1. — Spectroscopic redshift distributions to Kb = 21 for the various samples considered here. The BX/BM and BzK /SF distributions 
(hashed histograms) include sources from our sample in the GOODS-North field and overlap almost completely. The DRGs h ave a higher 
mean redshift of (z) = 2.48 ± 0.38 from our sample of J — Kg > 2.3 sources with 2: > 1 in all four fields of the optical survey iSteidel et alj 
120041) where we have complementary J and X-band imaging. The redshift distribution of DRGs with i n our sample (all of whic h are 
selected with the BX/BM or 2 3 LBG criteria) is similar to that found bv Ivan Dokkum etail <200l) . Ivan Dokkum et al.' C2003I), and 
^^rstcr Sclircibcr ct al. |2004j). The solid histogram shows the redshift distribution for BzK/PE galaxies from Daddi ct al. 1 2004a) and 
IDaddi et al... L2005.') . scaled down by a factor of 3 for clarity. The background shaded regions show the arbitrarily normalized redshift 
distributions for optically-selected BX and BM galaxies, and LBGs. 



extremely red colors which raay indicate passive evolu- 
tion. To address the former is sue, we have emplo yed 
the ''BzK'' selection criteria of IDaddi et al.1 lj2004a(l to 
cull objects in the GOODS-N field and directly compare 
with those selec ted on the observed optical properties of 
z > 2 galaxies. IDaddi et al.1 l)2004a|) define the quantity 
''BzK" : 

BzK ={z-K)-{B~ z); (3) 

star-forming galaxies with z > 1.4 are targeted by the 
following criterion: 

BzK > -0.2, (4) 

in AB magnitudes. Of the 1185 sources with > 3 a B, z, 
and K detections and K^ < 21, 221 satisfy Equation^ 
The surface density of BzK galaxies with Kg < 21 is 



~ 3 arcmin"^, similar to the surface density of BX/BM 
galaxies to a similar isTg-band depth. These star-forming 
BzK galaxies will be referred to as "BzK /SF" galax- 
ies, and their spectroscopic redshift distribution from our 
spectroscopic sample is shown in Figure^ Our deep near- 
IR imaging allows us to determine the redshift distribu- 
tion for BzK /SF galaxies with Kg > 20 (and which also 
satisfy the BX/BM criteria), and the results are shown in 
FigurelU The mean redshifts of the K^ < 20 and K^ > 20 
distributions are (z) = 2.13 ±0.22 and (z) = 2.03 ±0.41, 
respectively, and agree within the uncertainty. We note, 
however, that the BzK /SF criteria select Kg > 20 ob- 
jects over a broader range in redshift (1.0 < z < 3.2) than 
Ks < 20 objects. This reflects the larger range in BzK 
colors of Ks > 20 BzK /SF galaxies compared with those 
having Kg < 20. Additionally, the photometric scatter 
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in colors is expected to increase for fainter objects, so a 
broadening of the redshift distribution for BzK/SY ob- 
jects with fainter magnitudes is not surprising. 
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K<=20 BzK/SF ((z)= 2.13±0.22) 
K>20 BzK/SF ((z)= 2.03±0.41) 




Fig. 2. — Arbitrarily normalized spectroscopic redshift distribu- 
tion of BzK/SF galaxies in our spectroscopic sample to Kb ~ 22.5, 
with separate emphasis on Ks < 20 and Ks > 20 BzK/SF objects. 
The BzK/SF criteria select Kb > 20 objects over a broader range 
in redshift than Kb < 20 objects. 

We emphasize that we only know the redshifts for 
BzK/SF galaxies which also happen to fall in the 
BX/BM sample. In general, the true redshift distri- 
bution, No^^^^^ [z), of the BzK/SF sample will be 
broader than the distributions shown in Figure Q and 

Figure call them Nc'^^^^^{z)), which are effectively 
convolved with the BX/BM selection function. For ex- 
ample, the rapid dropoff in Nc^^^^^ [z) for z > 2.6 (Fig- 
ure 121) may simply reflect the dropoff in the BX selec- 
tion function for z > 2.6. However, the Nc^^^^^(z) 
we derive here is similar to that of the photometric 
redshift distribution of K20 galaxies from Daddi et aTl 
l)2004a(l . which is subject to its own systematic errors, 
suggesting that a reasonable approximation is to take 
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Near-IR Selection of Passively Evolving Galaxies 

In addition to the criteria above, several methods have 
been developed to select passively evolving high redshift 
galaxies by exploiting the presence of absorption or con- 
tinuum breaks in the SEDs of galaxies with dominant old 
stellar populations. The BzK selection criteria 

BzK <-Q.2 

z-K>2.5 (5) 

are designed to select passively evolving ga laxies at z > 
1.4 ijCimatti et al.l2004tlD'addi et al.l2004aj) . One galaxy 
which has a secure B-band detection, and an additional 
16 with _B-band limits, satisfy these criteria, implying a 
surface density of BzK/PF galaxies of 0.24 arcmin^^ to 
Kg = 21. Galaxies selected by their BzK colors to be 




3000 6000 9000 12000 
Rest-Frame Wavelength (A) 

Fig. 3. — Relative transmission of the z, J, Kb, and IRAC 3.6 /im 
filters at rest-frame wavelengths for z = 2.2. Also shown are typ- 
ical (unreddened) galaxy SEDs assuming constant star formation 
(CSF) and instantaneous star formation (r = 10 Myr) aged to 
1 Gyr. For redshifts z ~ 1.88 — 2.38, the J-band brackets the 
prominent Balmer and 4000 A breaks. 
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BX/BM 
BzK/SF ^ 
BzK/PE 
DRG ^3 



T 




T 



T 



Fig. 4. — Optical magnitude distributions for photometrically 
selected Kb < 21 galaxies in the BX/BM, BzK, and DRG samples. 
The solid vertical line denotes the magnitude limit for galaxies in 
the optically-selected (BX/BM) sample. Approximately 47% of 
DRGs (34/73) have TZ > 27.0 and are not shown in the figure. The 
distribution of BzK /PE galaxies has been arbitrarily normalized 
for clarity. 



passively-evolving are referred to as "BzK /PE" objects. 
The redshift distribution of BzK /PE galaxie s , taken 
fro m the spectroscopic samples of 'Dadd i et al.l lj2004a(l 
and lDaddi et al.l l)2005f) . shows that they mostly lie be- 
tween redshifts 1.4 ^ z < 2 (Figure^. We note that 
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Table 2. Sample Properties 



(-f'2.0-10kcv) (SFRx) 



Ka Range 


Sample 










(xlO**^ ergs s-i) 


(M0 yr-i) 


18.0 <Ks< 20.0 


BX/BM° 


11 


4 


10 (7) 


1.75 (1.80) 


7.13 ±0.88 (4.95 ±1.15) 


143 (99) 




BzK/SF 


77 


24 


45 (42) 


1.97 (2.01) 


5.39 ±0.46 (4.68 ±0.48) 


108 (94) 






14 


3 


11 


1.74 


Q (17 -L n 04 


61 




DRG 


20 


5 


14 


2.48 


5.26 ± 1.28 


105 


20.0 < -R's < 20.5 


BX/BM ° 


21 





21 


2.03 


4 89 + 79 


98 




B^ii'/SF 


57 


5 


56 (55) 


2.01 (2.01) 


3.96 ±0.40 (3.87 ±0.41) 


79 (77) 




















DRG 


20 


7 


13 


2.48 


2.32 ± 1.30 


46 


20.5 < Ks< 21.0 


BX/BM = 


27 





27 


1.99 


2.11 ± 0.56 


42 




rSzJx / or 


/ 


3 


82 (81) 


2.01 (2.01) 


2.32 ± O.oO (2.20 ± 0.30; 


4b (44 ) 




B^/C/PE 


3 


1 


2 


1.74 


2.07 ± 1.54 


41 




DRG 


33 


7 


26 


2.48 


1.90 ±0.86 


38 


21.0 <Ks< 21.5 


BX/BM <= 


31 





31 


2.04 


2.95 ± 1.01 


59 




BzK /SF 


99 


2 


97 


2.01 


2.76 ±0.28 


55 


21.5 <Ks< 22.0 


BX/BM = 


26 





26 


2.22 


1.43 ±0.77 


29 




Bzi^/SF 


148 





148 


2.01 


2.78 ±0.73 


56 


22.0 <Ks< 22.5 


BX/BM f 


93 





93 


2.04 


1.13 ±0.25 


23 




Bzit/SF 


77 





77 


2.01 


0.49 ±0.29 


10 



Note. — Values in parentheses are when we exclude all directly-detected X-ray sources, including ones which may be star-forming galaxies 
(Table EJ. 
''Total number of sources in sample. 

''Number of direct X-ray detections, corresponding to a minimum 3 a flux of /o. 5-2.0 kcV ^ 2.5 X 10^^^ erg s^^ cm^^. 
Number of stacked sources. 

''Mean redshift of stacked sample. Sources without a spectroscopic redshift are assigned the mean redshift of the sample to which they 
belong, where the mean redshifts for the sample are specified in Figure 111 

"^We only consider BX/BM galaxies which are spectroscopically confirmed to lie at redshifts z > 1, ^ 25% of the photometric sample of 
BX/BM galaxies. 

'Photometric BX/BM galaxies. 



we may be incomplete for the BzK/PE objects despite 
the very deep i?-band data considered here and these 
missing objects may be more easily selected using the 
J — Kg > 2.3 criteria discussed below (see also § 4.2.3). 

The J — Ks color probes the age-sensitive Balmer and 
4000 A breaks for galaxies with redshifts 2.0 ^ z < 4.5 
(Figure ISJ- The criterion 

J-i^s>2.3 (6) 



l|Franx et al.l l2003h can be used to select both passively 
evolving and heavily reddened star-forming galaxies with 
E{B—V) > 0.3. Galaxies satisfying this crit erion are also 
referr ed to as Dist ant Red Gala xies (DRGs; iFranx et al.l 
120031 van Dokku m et alJl2004|) . There are 62 galaxies 
with J — Kg > 2.3 that are detected in J, and an ad- 
ditional 11 have J-band limits. The observed surface 
density of DRGs is 1.01 ± 0.12 arcmin'^ to Kg = 21, in 
very good agreeme n t with the s urface density found by 
|van Dokk um et alJ l)2004f) and iForster Schreiber et al.l 
l|2004) . The spectroscopic redshift distribution of DRGs 
from the four fields of the optical survey where we have 
deep J and i^g-band imaging is shown in Figure ^ 
and is consistent with the redshift distributions found 



bv lvan Dokkum et alJ ll2004l) . [van Dokkum et al. (2003), 



yfi^ T -1.- -Y'r*»^ 

and IF5rster Schreiber et alJ (|200J). Star- forming and 
passively-evolving DRGs are referred to as "DRG/SF" 
and "DRG/PE", respectively. The depth of our J-band 
image implies that our sample of DRGs will be incom- 
plete for those with Kg > 21. Therefore, we have lim- 



ited ourselves to galaxies with Kg < 21 when comparing 
DRGs with BzK and/or BX/BM selected galaxies. We 
reconsider BX/BM and BzK galaxies with Kg > 21 as 
noted below. 

The optical magnitude distributions for galaxies with 
Kg < 21 are shown in Figure H The catalog of BX/BM 
galaxies is restricted to 7?, < 25.5. However, our optical 
imaging is significantly deeper {TZ = 27.5; 3 a), allowing 
us to extract optical magnitudes for galaxies much fainter 
than those in the BX/BM catalog. Most of those galaxies 
with TZ > 25.5 are DRGs. The nature of optically-faint 
DRGs is discussed in § 14.21 

For most of the analysis that follows, we either use only 
the spectroscopically confirmed sub-sample of BX/BM 
galaxies, or we apply our knowledge of the contamina- 
tion fraction of the photometric sample to deduce any 
inferred quantities. The small available spectroscopic 
samples using the near-IR criteria prevent us from ap- 
plying similar corrections when deducing properties for 
the near-IR samples. 

2.3. X-ray Data and Stacking Method 

One focus of this paper is to draw comparisons between 
galaxies selected by the techniques described above by 
using their stacked X-ray emission as a proxy for their 
bolometric SFRs. X-ray stacking allows us to deter- 
mine instantaneous bolometric SFRs in a manner that 
is independent of extinction and the degeneracies as- 
sociated with stellar population modeling. For exam- 
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pie, the average reddening of rest-frame UV selected 
galaxies of E{B — V) ^ 0.15 implies a column density 
of Nn ^ 7.5 X 10^° cm~^, as suming the Galactic cali- 
bration ijDiplasfc Savagelll994|) . Absorption in the rest- 
frame 2—10 keV band is negligible for these column 
densities. The X-ray data are taken fr om the Chandra 
2 Ms survey of the GOODS-N field llAlexander et all 
|2003|) . We made use primarily of the soft band (SB; 
0.5 — 2.0 keV) data for our analysis, but we also in- 
clude hard band (HB; 2.0 — 8.0 keV) data to examine 
the nature of directly-detected X-ray sources. The data 
are corrected for vignetting, exposure time, and instru- 
mental sensitivity in producing the final mosaicked im- 
age. The final product has an SB on-axis sensitivity of 
^ 2.5 X 10"^^ erg s~^ cm~^ (3 a), sufhcient to directly 
detect L2-10 kcv > 9.3 x 10^-^ ergs s~^ objects at 2 ~ 2, 
corresponding to an SFR of ^ 190 Mq yr^^. 

The stacking proced ure followed here is the same as 
that discussed in R.ed dv fc Steidell 1)20041) . Apertures 
used to extract X-ray fluxes had radii set to 2'.'5 for 
sources within 6' of the average Chandra pointing origin, 
and set to the 50% encircled energy ra dius for sources 
with off-axis angles greater than 6' (|Feigelson et al.l 
l2002t) . X-ray fluxes were computed by adding the counts 
within apertures randomly dithered by 0'.'5 at the galaxy 
positions. Background estimates were computed by ran- 
domly placing the same sized apertures within 5" of the 
galaxy positions, careful to prohibit the placing of a back- 
ground aperture on a known X-ray detection. This pro- 
cedure of placing random apertures was repeated 1000 
times. The mean X-ray flux of a galaxy is taken to be 
the average of all the flux measurements from the 0'.'5 
dithered apertures and the background noise is taken to 
be the dispersion in fluxes measured from the background 
apertures. We applied aperture corrections to the fluxes 
and assumed count rate to flux conversions based on th e 
results compiled in Table 7 of Alex ander et alJ l)2003D . 
a photon index T = 2.0, and a Galac tic absorption co l- 
umn density of A^h = 1.6 x lO^^^ cm^^ ||Stark et al.ll992|) . 
Poisson errors dominate the uncertainties in flux. 

3. RESULTS 

3.1. Direct X-ray Detections 

Of the 221 BzK/SF candidates with Ks < 21, 
32 (14%) have an X-ray counterpart within 1'.'5 
((Alexande r et al.ll200 3). with a 0.22% probabihty for 
chance superposition. The X-ray detection fractions are 
24%, 6%, and 26%, for the BzK/FE, BX/BM, and 
DRG samples, respectively, and are summarized in Ta- 
ble 12 Eleven of the 36 directly-detected BzK sources 
(32 in the BzK /SF sample and 4 in the BzK /FE sam- 
ple) would not have been detected with the sensitivity 
of the shall ower 1 Msec data in the GOODS-South fleld 
studied bv iDaddi et alJ l|2004a|) . After taking into ac- 
count the sensitivity difference, we find a direct X-ray 
detection rate comparable to that of Daddi et al. ( 2004a) 
of ^ 11%. Figure |5l shows the X-ray/optical flux ra- 
tios (log/x/Zn) for sources in all four samples (BX/BM, 
BzK/SF, BzK /FE, and DRG) directly-detected in the 
Chandra hard band (2 — 8 keV). Direct hard band detec- 
tions must be AGN if they are at 2; ~ 2, since starburst 
galaxies with no accretion activity are expected to have 
little flux at rest- frame energies of 6— 20 keV. Indeed, the 
X-ray/optical flux ratios for directly detected hard band 




10-16 10-'^ IQ-i" 10-'3 

Fj-Bkev (erg cm-2 s-') 

Fig. 5. — Optical/ X-ray flux ratios, defined as log/x//R = 
log/x + 5.50 + H/2.5 IHornschemeier et alJ l2001^. for all directly 
detected hard band X-ray sources in the BX/BM, BzK (SF/PE), 
and DRG samples. The abscissa is the observed hard-band flux, 
corresponding to rest-frame energies of 6 — 24 keV, and the ordi- 
nate is the observed Cousins R magnitude from the compilation 
of Barecr et al. (2003). Large circles denote sources with spec- 
troscopic rcdshifts 2 > 1, and are likely AGNs given their direct 
hard band detections. Sources with hard band detections but no 
redshift identification are shown by the small squares. Almost all 
sources have optical/X-ray fiux ratios between —1 < log /x//r < 1 
(shaded region), values commonly found for AGNs. Those with 
log fx/ fa < — 1 are confirmed interlopers, shown by the large open 
squares. Starburst galaxies generally have log/x//R < —1. 



sources lie in the region typically populated by AGNs 
(shaded area of Figure |3J). A smaller fraction of galaxies 
with direct hard band detections (and the three sources 
with the smallest log/x/Jn) are spectroscopically con- 
flrmed interlopers at z < 1. From Table it is easy to 
see that much of the AGN contamination in star-forming 
samples of galaxies (e.g, BX/BM, BzK/SF) occurs for 
magnitudes Ks < 20. 

Figure only shows those X-ray sources with hard 
band detections. Eleven additional sources had direct 
soft band (0.5 — 2.0 keV) detections, but no hard band 
detections. Of these 11, 5 sources have R > 22.0, 
/0.5-2.0 kov < 0.1 X 10-15 erg cm'^ s-\ and log/x//R < 
— 1, indicating they may be starburst galaxies. These 5 
sources and their properties are summarized in Table |3| 
T hree of th e five sources have spectra taken by us or 
bv lBarger et al.l l)20fK^ indicating no obvious AGN spec- 
tral features. These sources may be rapidly star-forming 
galaxies and for fairness we include them in the stacking 
analysis as indicated below and in Table 13 

It is interesting to also consider the rest-frame near-IR 
properties of the directly detected X-ray sources as indi- 
cated by their Spitzer/IRAG colors. Figure shows the 
3.6 — 5.8 fj,m color as a function of 3.6 fim magnitude 
for all samples considered here. There is a clear seg- 
regation in the IRAC colors of X-ray detections where 
they show, on average, brighter IRAC magnitudes and 
redder IRAC colors when compared with the colors of 
star- forming galaxies in the BX/BM and BzK/SF sam- 
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Table 3. Possible Star- Forming Direct X-ray Detections 



(J2000.0) (J2000.0) 



Sample 



fo.5 



2.0 keV 
erg cm~ 



(Vega mag) 



(mag) 



log/xZ/R*' 



12:36:21.95 


62:14:15.5 


BzK/SF 


1.38 


0.02 


18.95 


23.1 


-1.96 


12:36:52.75 


62:13:54.8 


BX/BM 


1.36 


0.03 


19.54 


22.1 


-2.18 


12:36:53.46 


62:11:40.0 


BX/BM;_Bzi<'/SF 




0.11 


18.65 


22.7 


-1.38 


12:36:56.89 


62:11:12.1 


BzK/SF 




0.02 


20.50 


23.8 


-1.68 


12:37:03.70 


62:11:22.6 


BX/BM;BzK/SF 


1.72 


0.04 


19.92 


23.4 


-1.54 



''Soft-band fluxes are from I Alexander et alj 120031) and Cousins R magnitudes are from lBarger et all 1200 
''Defined as log/x//R = log/x -I- 5.50 + i?/2.5 IHornschemeier et alJl2001^ . 
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Fig. 6. — Spitzer/TRAC 3.6 — 5.8 ^tm color versus 3.6 fim magni- 
tude (in AB units) for all samples considered here, with emphasis 
on directly detected X-ray sources (large circles). These direct 
detections generally have brighter IRAC magnitudes and redder 
colors than star forming BX/BM and BzK/SF galaxies, likely due 
to thermal continuum from circumnuclear dust proximate to the 
AGN. 



pies. Such a trend might be expected if the rest-frame 
near-IR hght from the X-ray sources is dominated by 
thermal continuum from circumnuclear dust heated by 
the AGN. The increase in flux density across the IRAC 
bands for AGN has been seen for ERO samples at red- 
shifts z ~ 1 — 3 ijFraver et al.l200^ , similar to what is ob- 
served here. Finally, the 5 objects listed in Tabled have 



'7^3. 6^111-5. 8^111 



-0.35 to 0.35 and ma.gum = 20.2 - 22.0 



lying in the same region of IRAC color space as some of 
the star-forming BX/BM and BzK/SF candidates. 

At times in the following analysis, we also consider 
submillimeter galaxies and their relation to optical and 
near-IR selected objects. These heavily star-forming ob- 
jects are generally associated with directly-detected X- 
ray sources, and we reconsider the X-ray emission from 
these sources as pointed out below. Unless otherwise 
stated, however, we have excluded all directly-detected 
hard band X-ray sources from the analysis under the as- 
sumption that their X-ray emission is contaminated by 
AGN. 
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0.04 



' ' ' I ' ' ' ' I ' ' ' ' I 
BX/BM: ((z-K)^b)=1.07±0.69 
BzK/SF: ((z-K)^B>=1.61±0.61 
.BzK/PE: ((z-K)jb) = 2.80±0.26 
DRG: ((z-K)^b) = 2-49±0.78 




Fig. 7. — z — K color distribution for BX/BM, star-forming 
BzK, and DRG galaxies to Kb = 21. The mean z — K color of 
galaxies becomes redder for the BX/BM to BzK to DRG samples. 
I3X/BM selection is more efficient in selecting objects with blue 
{z — K)ab ^ 1 and DRG selection is more efficient in selecting 
objects with very red (z — K)\b > 3. The BzK criteria spans 
the middle range of (z — K)\^ color. The small solid histogram 
shows the arbitrarily normalized distribution in (z — -R')ab color 
for passively-evolving BzK galaxies. 



3.2. Overlap Between Samples 

Galaxies selected solely by the presence of some un- 
obscured star- formation (BX/BM selection), and those 
selected by some combination of stellar mass and star- 
formation (DRG and BzK selection) can be distin- 
guished by their observed near-IR color distributions 
(FigureEI). The mean {z-K)ab color for BX/BM galax- 
ies is ^ 0.54 mag bluer than the BzK/SF sample, just 
within the 1 a dispersion of both samples. This differ- 
ence in average (z — K)ab color between BX/BM and 
BzK/SF galaxies partly stems from the fact that the 
width of the BzK selection window below (z — K) ab = 1 
narrows to the point where photometric scatter becomes 
increasingly important in determining whether a galaxy 
with blue colors (i.e., (z — K)ab < 1) is selectable with 



Optical/Near-IR Selected Galaxies at Z ~ 2 



9 



the EzK/SY criteria^. On the other hand, the BX/BM 
criteria are less efficient than BzK selection for galaxies 
with {z — K)as, ^ 1-6. BX/BM galaxies with red near-IR 
colors are systematically fainter in the optical than those 
with blue near-IR colors (Figure |H1), reflecting both the 
correlation between TZ and z as these filters lie close in 
wavelength, as well as the Kg < 21 limit adopted in Fig- 
ure [S] Therefore, the optical catalog limit oiTZ = 25.5 
would appear to exclude from the BX/BM sample those 
galaxies with {z — K)p^-q > 3 (Figure |H1)- As we show 
§ ^ the exclusion of (z — K)ab ^ 3 galaxies by opti- 
cal selection is not a fault of the criteria themselves: the 
TZ = 25.5 limit is imposed so that spectroscopic followup 
is feasible on the candidate galaxies. Rather, the exclu- 
sion of (z — K)ab ^ 3 galaxies from optical surveys sim- 
ply reflects a fundamental change in the star formation 
properties of such red galaxies. 
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Fig. 8. — {z ~ K)\B versus TZ for photometric BX/BM galaxies 
(squares) to Kb = 21, showing that optically selected galaxies with 
red near-IR, colors are optically fainter on average than those with 
blue near-IR colors. This effect is due to the Kb = 21 limit as 
well as the correlation between TZ. and z-band magnitude as the 
two filters lie close in wavelength. Objects with {z — K)ab 
2.6 are missed by BX/BM selection as they fall below the TZ = 
25.5 BX/BM catalog limit (dashed vertical line). Also shown are 
BX/BM sources with direct X-ray detections (large circles) and 
spectroscopically confirmed interlopers (triangles). 

Separately, DRGs have a very red ((z — K)ab) = 
2.49 ± 0.78. Approximately 10% of z > 1.4 BzK /SF 
galaxies also have J — Ks > 2.3 (Figure El), similar to 
that found bv lDaddi et al.l l)2004a|) . The fact that there 
is some, albeit small, overlap between the BzK /SF and 
DRG samples is not surprising since the two criteria can 
be used to target reddened galaxies and both have red- 
shift distributions that overlap in the range 2.0 < z < 2.6 
(Figure The DRG fraction among BzK selected 
galaxies does not change appreciably if we add in the 

None of the selection criteria considered here have boxcar se- 
lection functions in either color or redshift space due to various 
effects, including photometric errors. The effect of this is to sup- 
press the efficiency for selecting objects whose intrinsic colors lie 
close to the edges of the selection window. 
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Fig. 9. — J - Kb color versus Kb for BX/BM (filled squares), 
BzK (open squares), and DRG (circles) samples to Kb = 21. The 
hashed horizontal line denotes the J — Kb = 2.3 limit. DRGs with 
limits in J-band are indicated by the triangles. Approximately 
5% of DRGs satisfy the BX/BM criteria, but the fraction rises to 
~ 12% if we include those selected by the z ~ 3 LBG criteria. 




Fig. 10. — (a) Fraction of BzK /SF sources which are optically 
selected with the BX/BM criteria when including (squares) and 
excludin g (c ircles) directly detected X-ray sources that are likely 
AGN (§E}; (b) Fraction of photometric BXs and BMs that are 
BzK /SF selected (squares) and the fraction of BXs and BMs with 
confirmed redshifts z > 1.4 that are BzK /SF selected (circles). 



BzK/PE sources — only 5 of 17 BzK/PF galaxies have 
J—Ks > 2.3 — as the BzK/PF galaxies are mostly at red- 
shifts lower than the DRGs {z <2). Finally, we note that 
DRGs include objects with much redder (z — K) ab colors 
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Fig. 11.— BX/BM colors of BzK/SF galaxies (circles) and 
DRGs with Ti, < 25.5, excluding direct X-ray detections. Squares 
denote DRGs with > 5 cr detections in C/n, G, and 7?,; triangles 
denote DRGs with 5 u limits in (7n. AH points are for galaxies with 
< 21. Also indicated are the BX/BM selection criteria for z ^ 2 
BMs and BXs, as well as the z ^ 3 LBG criteria of Steidel et al.l 
12003) . Approximately 40% of our DRG sample galaxies have TZ < 
25.5. Of all DRGs, including those not shown in the figure and 
which have Ti > 25.5, ~ 12% can be selected using the BX, BM, 
or LBG criteria. Note the number of BzK and DRG galaxies that 
lie very close (e.g., within < 0.2 mag) of the BX/BM selection 
windows. 
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Fig. 12. — BzK colors of BX/BM galaxies (empty circles) and 
DRGs (filled squares) to Kb = 21. Large filled circles denote 
BX/BM objects with spectroscopically confirmed redshifts 2 < 1 
(interlopers), most of which fall outside the BzK/SF selection win- 
dow. Also shown is the expected region of color space for passively 
evolving 2 > 1.4 galaxies fDaddi et al. 2004a). Note the number 
of BX/BM galaxies that lie very close (e.g., within < 0.2 mag) of 
the BzK selection window. DRGs with _B-band limits, shown by 
the triangles, cluster in the region expected for passively evolving 
2 > 1.4 galaxies. 



than found among BX/BM and BzK/SF/PE galaxies, 
i.e. those with {z — K)ab > 3. The absence of these 
galaxies from star-forming selected samples is discussed 
in S IO 

We can directly quantify the overlap between BzK/SF 
and BX/BM galaxies. Figure ^] shows the fraction of 
BzK/SF galaxies satisfying the BX/BM criteria (left 
panel) and the fraction of BX/BM galaxies satisfying the 
BzK/SF criteria (right panel). Most of the contamina- 
tion of the BzK/SF sample (that we know of) is from 
X-ray detected AGN 13. 1|) . while most of the contami- 
nation of the BX/BM sample is from low redshift inter- 
lopers (Table . Both sources of contamination tend to 
occupy the bright end of the K-hand apparent magni- 
tude distribution. We also show the overlap fractions in 
Figure E| excluding X-ray detected AGN and interlop- 
ers. The BX/BM criteria recover an increasing fraction 
of BzK/SF selected sources proceeding from Kg < 20 
galaxies 60% recovery fraction) to Kg ~ 21 galaxies 
(^ 80% recovery fraction) after excluding directl y de - 
tected X-ray sources which are likely AGN (see § I3.1|> . 
Conversely, the BzK/SF criteria recover ~ 80% of spec- 
troscopically confirmed BX/BM galaxies at z > 1.4, 
and are evidently effective at recognizing most of the 
BX/BM low redshift interlopers that tend to occupy the 
bright end of the if -band apparent magnitude distribu- 
tion. This result stems from the fact that low redshift 
interlopers tend to have bluer colors than necessary to 
satisfy the BzK/SF criteria. 

Figures and ^] show that a significant portion of 



BzK/SF galaxies missed by BX/BM selection, and con- 
versely, have colors that place them within < 0.2 mag of 
the selection windows, which is comparable to the pho- 
tometric uncertainties. The BX/BM criteria likely miss 
some BzK/SF galaxies not because of some failure of 
the criteria, but because we cannot measure photometry 
with infinite precision. The trend from lower (60%) to 
higher (80%) recovery rate shown in Figure llOh . reflects 
the fact that a greater percentage of Kg < 20 BzK/SF 
galaxies have redder G — TZ colors (when compared with 
Kg > 20 BzK/SF galaxies) than required to satisfy the 
BX/BM criteria (Figure [TT|l . There are some BzK/SF 
galaxies which have very red G — TZ > 0.8 colors. As 
we show in § 21 these red G — TZ galaxies would have 
an average bolometric SFR similar to BzK/SF galaxies 
with bluer G — TZ colors if they are at similar redshifts, 
z ^ 2. Therefore, if these red objects are at z ^ 2, then 
the correlation bet ween G — TZ an d reddening, as quan- 
tified by the ICalzetti et al.l l|20fl(1l) law, would appear to 
fail. Photometric scatter will also reduce the effective- 
ness of the BzK criteria in selecting BX/BM galaxies 
(Figure El- We can account for most of the photomet- 
ric in^om^l^tenessjasing the more sophisticated analysis 
of lReddv et alJ ll200l . 

Our deep K-hand data allow us to investigate the ef- 
ficiency of BzK/SF selection to fainter K magnitudes 
than previously possible. Figure E| shows the BzK 
colors of BX/BM galaxies with spectroscopic redshifts 
1.4 < z < 2.6 for three bins in Kg magnitude. The 
BzK/SF criteria were designed to select relatively mas- 
sive galaxies with Kg < 20, but they become slightly 
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less efficient in culling Kg > 21 galaxies: 10 of 49 
(~ 20%) BX/BM galaxies with spectroscopic redshifts 
1.4 < z < 2.6 and > 21 do not satisfy the BzK/SF 
criteria. Furthermore, we note that ~ 11% (61/544) of 
BX/BM candidates that fall in the region with 7^-band 
data are undetected to Kg = 22.5 (3 a). The iiT-band 
limits for these galaxies suggests they are younger star- 
forming systems with (z — K)ab ^ 1, below which the 
BzK /SF criteria drop in efficiency, as discussed above. 
We remind the reader that many of the _BzA'/SF objects 
not appearing in the BX/BM sample may be missed by 
the BX/BM criteria simply because of photometric er- 
rors. BX/BM galaxies missed by the BzK/SF criteria 
may be missed not because of intrinsic differences in the 
objects, but simply because of photometric scatter or 
because of the difficulty in obtaining very deep i^g-band 
data. 
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Fig. 13. — BzK colors of spectroscopically confirmed BX/BM 
galaxies with redshifts 1.4 < z < 2.6 for different bins in Kb mag- 
nitude. BzK/SF selection appears to miss an increasing fraction 
of Kb > 21 galaxies in this redshift range due to the narrowing 
BzK/SF selection window for objects with bluer (z — K)\^ col- 
ors. In addition, ~ 19% of BX/BM candidates have Kb > 22.0 and 
are not shown. 

Turning to J — Kg > 2.3 galaxies, we show the opti- 
cal colors of DRGs with TZ < 25.5 in Figure ITTl and the 
near-IR colors of the 74% of DRGs with z-band detec- 
tions in Figure IT^ The optical criteria are particularly 
inefficient in selecting J— Kg > 2.3 sources: 9 of 73 (12%) 
of DRGs in the GOODS-N field satisfy BX, BM, or LBG 
selection. This fra ction is similar t o the overall detec- 
tion rate found bv lErb et al.l l)2fl05D for the 4 fields in 
the z ~ 2 optical survey with deep J and if-band data. 
The LBG criteria can be used to select some DRGs since 
z ~ 2 galaxies with a Calzetti reddened constant star 
formation SED with E{B ~ V) > 0.3 are expected to lie 
in the color space occupied by z ~ 3 LBGs. A greater 
fraction (- 30%) of DRGs satisfy the BzK/SF criteria 
(Figure [T^ since these criteria select objects with devel- 
oped spectral breaks and with redshifts that fall in the 
range probed by DRG selection (see Figure^. We note 



that DRGs with i3-band limits cluster in the region of 
color space expected for passively evolving z > 1.4 galax- 
ies (Figure El . These DRGs have little, if any, current 
star formation (see § I4.2|l . 

3.3. Stacked X-ray Results 

The X-ray data are not sufficiently sensitive to de- 
tect individual galaxies with SFR< 190 Mq yr~i (3 a). 
We can, however, stack the X-ray data for subsets of 
galaxies below the sensitivity threshold to determine 
their average X-ray emission. The influence of AGN 
in any X-ray stacking analysis is a concern. The soft- 
ness of a stacked signal provides some circumstantial ev- 
idence for X-ray emission due primarily to star forma- 
tion (e.g., f van Dokkum ct al. 2004, ' Daddi et a i''2004al 
iLaird et al.i „2005 ). UV line signatures and radio emis- 
sion can provide additional constraints on the presence 
of AGN fe.g.. IR,eddv fc Steide]ll200l . We typically re- 
moved all directly-detected X-ray sources from the op- 
tical and near-IR samples before running the stacking 
simulations, except as noted below and in Table when 
considering X-ray detected sources which may be star- 
forming galaxies. Our method of excluding other X-ray 
detected sources ensures that luminous AGN do not con- 
taminate the stacked signal. Indirect evidence suggests 
that less luminous AGN do not contribute significantly 
to the stacked signal. First, the stacked signal has no 
hard band (HB; 2 — 8 keV) detection indicating that the 
signal is softer than one would expect with a significant 
AGN contribution. Second, the availability of rest-frame 
UV spectra for many of the BX/BM objects provides 
an independent means of identifying AGN. There is one 
source whose spectrum shows high-ionization emission 
lines in the rest-frame UV, but no X-ray detection in the 
Chandra 2 Ms data. Removing this X-ray faint AGN 
source does not appreciably af f ect th e stacked X-ray flux. 
In addition, iReddv fc Steidell l|2004D examined the very 
same BX/BM dataset used here and found a very good 
agreement between dust-corrected UV, radio, and X-ray 
inferred SFRs for the sample, suggesting star formation 
as the dominant mechanism in producing the observed 
multi- wavelength emission. Finally, the local hosts of low 
luminosity AGN have stellar populatio ns characteristic of 
passi vely ev olvin g early-type galaxies ijKauffmann et afl 
we show that passively evolving galax- 
ies at z ^ 2 have little or no detectable X-ray emission, 
implying that low level accretion activity in these sys- 
tems does little to alter the X-ray emission relative to 
that produced from star formation. The absence of X- 
ray emission from these passively evolving galaxies also 
suggests that low mass X-ray binaries contribute little 
X-ray emission in star-forming galaxies when compared 
with the emission produced from more direct tracers of 
the current star formation rate, such as high mass X-ray 
binaries. 

Stacking results for the samples (to Kg — 22.5) are 
summarized in Figure 174] and Table |21 The left panel of 
Figure ITl] includes all photometrically-selected BzK/SF 
and DRG galaxies, and all spectroscopically-conffimed 
z > 1 BX/BM galaxies. The right panel includes 
only those BzK/SF galaxies with spectroscopic redshifts 
z > 1, all of which also satisfy the BX/BM criteria. All 
direct X-ray detections have been excluded in making 
Figure [T^ The distributions do not change appreciably 
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if we only consider the X-ray flux of BzK /SY galaxies 
spectroscopically confirmed to lie at z > 1 (Figure ^^). 
Removing the one spectroscopic z > 1 AGN undetected 
in X-rays does little to change the X-ray luminosity dis- 
tributions. The luminosity distributions agree well be- 
tween the three samples over a large range in Ks magni- 
tude, with Ks < 20 galaxies exhibiting the largest X-ray 
luminosities by a factor of two to three when compared 
with fainter > 20.5 galaxies. 
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Fig. 14. — Stacked X-ray luminosity versus Kb magnitude for 
BX/BM galaxies with 2 > 1 (circles), BzK /SF galaxies (triangles), 
and DRGs (squares), excluding all directly detected X-ray sources. 
The left and right panels show the distributions for photometric 
and spectroscopically confirmed z > I BzK/SF galaxies, respec- 
tively. Sources without a spectroscopic redshift were assigned the 
mean redshift of the sample which they belong to, according to 
Figure In all cases, we find the distributions consistent within 
the errors. We also find that Kb < 20 galaxies have average X-ray 
luminosities that are a factor of 2 — 3 times higher than that of 
Kb > 20.5 galaxies. 



4. DISCUSSION 

In this section we first present the X-ray inferred 
average bolometric SFRs for galaxies in the BX/BM, 
BzK /SF /PE, and DRG samples, and compare our re- 
sults with other X-ray stacking analyses. Unless stated 
otherwise, we exclude hard band X-ray AGN sources 
from the analysis of the SFRs. The SFRs are interpreted 
for galaxies as a function of their near-IR colors and we 
assess the ability of optical surveys to single out both 
heavily reddened and massive galaxies. We identify pas- 
sively evolving galaxies at z ^ 2 from their red near-IR 
colors and discuss plausible star formation histories for 
these galaxies using the X-ray data as an additional con- 
straint. Finally, we discuss the contribution of BX/BM, 
BzK/SF, and DRG galaxies to the star formation rate 
density at z '--^ 2, taking into account the overlap between 
the samples and their respective redshift distributions. 

4.1. Star Formation Rate Distributions 



Star Formation Rates and Comparison with Other 
Studies 

We esti mated the SFRs for galaxies in our samples 
using the 'Ranalli e t al.l l|2003fl calibration between X- 
ray and FIR luminosity. This calibration reproduced 
the SFRs based on independent star forma tion tracers 
for z - 2 galaxies ijReddv fc Steidell l2004|) . so we are 
confident in using it here. The SFR distributions for 
BX/BM, BzK, and DRG galaxies are shown in Fig- 
ure ^1 where we have added the 5 directly-detected soft 
band X-ray sources in Table |2| that may be star- forming 
galaxies. The SFRs are summarized in Table El The 
mean SFR of K^ < 20 galaxies is - 90 - 140 Mq yr'^ 
and is a factor of 2 — 3 tim es larger than galaxi es with 
Kg > 20.5. For comparison, iDaddi et all t2QQ4af ) found 
an average SFR of K20 galaxies in the GOODS-South 
field of 190 Mq yr~^ (including one likely star-forming 
galaxy directly detected in X-rays). This is somewhat 
higher than our value of 110 M© yr^^ for A's < 20 
BzK/SF galaxies. This discrepancy could simply re- 
sult from field-to-field variations, small number statis- 
tics, or the lower sensitivity of the X-ray data in the 
GOODS-South field compared to GOODS-North. With 
the Chandra 2 Ms data, we are able to exclude directly- 
detected X-ray sources down to a factor of two lower 
threshold than was possible with the 1 Ms data in the 
GOODS-South field. If we add back those K^ < 20 X- 
ray BzK/SF galaxies that would have been undetected 
in the 1 Ms data to the stacking analysis, we obtain 
an average SFR of 160 Mq yr~^, more in line with the 
IDaddi et alJ l|2004aj) value of 19 M© yr-^ 

A similar stacking analysis bv iRubin et al.l 1)20041) indi- 
cates that Ks < 22 DRGs have SFRs of ~ 280 Mq yr^^, 
corrected for the difference in SFR calibration used here 
and in Rubin ct al. (2004). This very high value is likely 
a result of the shallow X-ray data (74 ksec) considered in 
that study; the depth of their X-ray data precludes the 
removal of most of the X-ray sources that are directly 
detected in the 2 Ms X-ray survey. If we include those 
X-ray sources that would have been undetected in the 
74 ksec data, assuming a mean redshift of (z) — 2.4, the 
average SFR for the DRGs with K^ < 21 is 250 M© yr"!. 
Therefore, much of the difference in the SFRs can be at- 
tributed to unidentified AGN in the shallower X-ray sur- 
veys contaminating estimates of the star formation rate. 
If Figure is any indication, then adding DRGs with 
21 < ATs < 22 to the stack would decrease this average 
SFR. Variance of the fraction of DRG/PE galaxies be- 
tween fields may also affect the average SFRs: a greater 
fraction of DRG/PEs in the GOODS-N field, - 25% 
(3/13) of which have K^ < 20, will lead to a lower av- 
erage SFR for Ks < 20 DRGs. As we show in § 4.4.2, 
there are clearly some number of very reddened galax- 
ies with large SFRs (e.g., submillimeter galaxies) among 
DRGs (and among the BX/BM and BzK/SF samples). 
Regardless, these calculations underscore the importance 
of factoring in the differing sensitivity limits of the vari- 
ous X-ray surveys before comparing results. The strong 
dependence of SFR with Ks magnitude (Figure I15|) also 
suggests that fair comparisons of the SFRs of galaxies 
selected in different surveys can only be made between 
objects with similar rest-frame optical luminosities. 

Our analysis is advantageous as we are able to compare 
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Fig. 15. — St ar formation rates i nferred from the X-ray lumi- 
nosity using the IRanaUi et al.l 120031) caUbration. We have added 
the 5 directly-detected soft band X-ray sources that may be star- 
forming (Table |3J to compute the average SFR. Circles, triangles 
and squares denote BX/BM, BzK, and DRG samples, respectively. 
The BX/BM points are for spectroscopically confirmed ^ > 1 
galaxies, and we have assumed the mean redshifts for the BzK /SF 
and DRG samples as in Figure^] Uncertainties in the star forma- 
tion rates are dominated by scatter in the X-ray/FIR relation and 
the dispersion in the stacked X-ray estimates. The dashed hori- 
zontal line denotes the average SFR for the en tire spectroscopic 
(z> 1) BX/BM sample of ~ 50 Mq yr"! from lReddv fc Steidell 
12004) . Galaxies with Kb < 20 have inferred star formation rates 
a factor of 2 — 3 higher than for Kg > 20.5 galaxies. 



the SFRs of galaxies within the same field, employing 
the same multi-wavelength data (to the same sensitiv- 
ity level) and the same photometric measurement tech- 
niques, for a consistent comparison. The inferred average 
SFRs of BzK /SF and DRG galaxies are remarkably sim- 
ilar to that of optically selected galaxies once the samples 
are restricted to similar Kg magnitudes. The previously 
noted discrepancies in X-ray inferred SFRs of BzK /SF, 
DRG, and BX/BM galaxies are therefore likely a result 
of a mismatch between X-ray survey limits and near- 
IR magnitude ranges. Field-to-field variations may also 
partly account for the previously observed discrepancies. 

Dependence of SFR on (z — K)ab Color 

We began our analysis by noting the differences be- 
tween the (z — K)ab color distributions of optical and 
near-IR selected galaxies (Figure [Tj). Figure ^1 indi- 
cates that despite these near-IR color differences, the 
BX/BM, BzK, and DRG galaxies have very similar av- 
erage SFR distributions as a function of Kg magnitude. 
Another pr oxy for stellar mass i s the (z—K) ab orTZ—Kg 
color fe.g.. iShaolev et ani200,'TD as it directly probes the 
strength of the Balmer and 4000 A breaks. Figure El 
shows the inferred average SFRs of optical and near-IR 
selected galaxies as a function of their (z — K)ab color, 
excluding all directly-detected X-ray sources^. Within 

^ Adding the 5 sources in Table U] does not appreciably affect 
Figure 



any single sample, objects with red {z—K)ab colors up to 
{z — K) AB ^ 3 have the largest SFRs. The red (z — K) ab 
color for these objects with high SFRs likely results from 
a developed spectral break (due to an older stellar pop- 
ulation) combined with the effects of dust. In fact. Fig- 
ure El illustrates the tendency for BX/BM objects with 
spectroscopic redshifts z > 1 and red {z — K) ab colors to 
have larger attenuation, as parameterized by E{B — V), 
on average, than those with bluer {z — K)ab colors. The 
turnover in the inferred SFR around {z — K)ab 3 is 
discussed in the next section. 
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Fig. 16. — Star formation rates of Kb < 21 galaxies as a func- 
tion of (z — K)a^ color. Symbols are the same as in Figure [TE\ 
The small solid triangles denote the average SFR inferred for BzK 
galaxies that are selected to be passively evolving (BzK /PE galax- 
ies) and the large inverted triangle indicates the limit in SFR found 
for the 13 DRGs with (z — K)ab > 3.0. Note the turnover in in- 
ferred SFR at (z — K)ab ~ 3. Also indicated are the average Kb 
magnitudes for sources in each bin of (z — K)a^ color. 



Figure El suggests that optically selected BX/BM 
galaxies may have systematically higher SFRs than BzK 
and DRG galaxies with similar {z — K)ab colors, per- 
haps indicating that the stacked sample for BzK/SF and 
DRG galaxies includes those which are passively evolv- 
ing. This may be particularly true of DRGs, quite a 
few of which only have B-band limits and which cluster 
in the BzK color space occupied by passively evolving 
galaxies (Figure I12|) . We can assess the dispersion in 
SFRs by separately stacking galaxies that are expected 
to be currently star-forming based on their colors and 
those that are not. For example, the average SFR com- 
puted for BzK/SF galaxies that are not selected by the 
BX/BM criteria (of which ~ 70% are within 0.2 mag of 
the BX/BM selection windows) is ^ 70 Mq yr~^, com- 
parable to the average SFR of all BzK/SF galaxies with 
Kg < 21. In summary, star-forming BzK galaxies have 
similar SFRs regardless of whether or not they satisfy 
the BX/BM criteria. The BX/BM criteria miss some 
BzK/SF galaxies either because of photometric scatter 
or because their optical colors are not indicative of their 
reddening. It is worth noting that photometric scatter 
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Fig. 17. — Attenuation, as parameterized by the rest-frame UV 
spectral slope, E{B — V), as a function of (z — K)ab color for 
spectroscopically confirmed BX/BM galaxies with z > 1. Errors 
in E(B — V) represent the 1 a dispersion of E(B — V) values for 
each bin of (z — K)a_b- 



works both ways: some sources with intrinsic colors sat- 
isfying the BX/BM criteria will be scattered out of the 
BX/BM selection windows and some whose colors do not 
satisfy the BX/BM criteria will be scattered into the 
BX/BM selecti on windows, althoug h the two effects may 
not equilibrate ijReddv et al.ll2005j) . The incompleteness 
of the BX/BM and BzK /SF criteria with respect to all 
star-forming galaxies at z ~ 2 simply reflects our inabil- 
ity to establish perfect selection criteria immune to the 
effects of photometric scatter and SED variations while 
at the same time efSc iently excluding interlopers (e.g., 
lAdelberger et al]l2004|) . However, the advantage of spec- 
troscopic optical surveys is that their s election functions 
can b e quantified r elatively easily fe.g.. lAdeiberger et ahl 
l2()()llR,eddv et 

Optical Selection of Reddened Star-Forming Galaxies 

Naively, one might interpret the inferred SFRs as a 
function of (z — i^)AB color combined with the results 
shown in Figure|Slto suggest that BX/BM selection may 
miss the most actively star-forming galaxies at z ^ 2. 
We can interpret the similarity in SFRs of BX/BM and 
BzK /SF galaxies in the context of their reddening, as 
parameteriz ed by their r est-fram e UV spectral slopes, 
E{B - V). iDaddi et all ^QOAa) show that the redden- 
ing vector is essentially parallel to the BzK limit de- 
fined by Equation^ implying that BzK selection should 
be sensitive to galaxies with higher extinction (and pre- 
sumably higher star formation rates) than found among 
BX/BM selected galaxies (i.e., E{B -V) > 0.3). How- 
ever, the similarity in the average SFRs of BX/BM and 
BzK /SF galaxies suggests several possibilities. First, 
we noted above that BzK /SF galaxies not selected by 
the BX/BM criteria have similar SFRs as those which 



do satisfy the BX/BM criteria*^. Consequently, BzK /SF 
galaxies with large SFRs that do not satisfy the BX/BM 
criteria because they truly have E{B — V) > 0.3 may 
not exist in sufficient numbers to significantly change 
the average SFRs for all BzK/ SF galaxies which do 
not satisfy the BX/BM criteria. lAdelberger fc Steidell 
()2000l) and lLaird etal] l|2005() find that optically selected 
galaxies with z > 1 show no correlation between their 
rest-frame UV luminosities and their obscuration, im- 
plying that on average the redder (more obscured) galax- 
ies have higher bolometric SFRs than galaxies with less 
reddening. Therefore, the similarity in the average X- 
ray inferred SFRs of BX/BM and BzK /SF galaxies sug- 
gests that there are not large numbers of galaxies with 
E{B — V) > 0.4 (i.e., if there were a large number of such 
heavily reddened objects, their bolometric SFRs would 
imply X-ray luminosities large enough to be directly de- 
tected in the soft band X-ray data, and very few likely 
star-forming galaxies at these redshifts are directly de- 
tected in the soft band). 

Second, studies of the UV emission from submillime- 
ter galaxies (SMGs) suggest that heavily reddened galax- 
ies may have similar rest-frame UV spectral properties, 
such as their range in E{B — V), as those which are 
forming stars at modest rates, implying that the cor- 
relation b etween E(B — V) a nd bolometric SFR (e.g.. 
from the iMeurer et all Il99fll and iCalzetti et all 1200(1 
laws) breaks d own for the most actively star-forming 
galaxies (e.g., iChanman et al. 2005). Table 01 sum- 
marizes the properties of the 9 known radio-selected 
SMGs with S'sso^m > 5 mJy in the GOODS-N field 
that overlap with our near-lR imaging l|ChaDman et alJ 
l200iTlWang et alJ2004D . Also Hsted are t wo (pair) sources 
wit h SnF,niim ^ 4 mJy t aken from Chapm lm et al.l l|200,'Tf) 



and lWang et all ll200l) . Of the 7 SMGs with redshifts 
1.4 < z < 2.6, 3 satisfy the BX/BM criteria. The de- 
tection rate of ~ 40% is similar to th e detection rate 
of SM Gs with BX/BM colors found bv iChaoman et alJ 
()2005j) . The mean bolometric luminosity of the 5 SMGs is 
(iboi) ~ 9 X 10^^ Lq as inferred from their submillimeter 
emission, correspond ing to an SFR of ^ 1500 M0 yr~^ 
using the iKennicutti il 1 99.S|i relation^. Despite their large 
bolometric luminosities, the three submillimeter galax- 
ies with redshifts in our sample have dust-corrected UV 
SFRs of 14—28 Mq yr~-^. In these cases, the UV emission 
may come from a relatively unobscured part of the galaxy 
or ma y be scattered ou t of the optically-thick dusty re- 
gions (jchanman et al.ll200,'iD . The BzK /SF criteria cull 
5 of the 7 SMGs with redshifts 1.4 < z < 2.6. There- 
fore, at least in the small sample of SMGs examined here 
(irrespective of their X-ray properties), the BX/BM and 
BzK /SF samples host an approximately equal number 
of SMGs. Finally, as we show below, galaxies with the 
most extreme (z — K)ab colors (i.e., (z — K)ab > 3) are 
red not because they are obscured by dust, but because 
they have little or no current star formation. It is there- 
fore not surprising that such objects are not identified by 

Those BzK /SF galaxies with G - 7?. > 1 and blue Un-G <1 
have optical colors that are similar to the c olors expec ted for lower 



redshift {z ^l) galaxies (e.g.,Elelbereer eOnEHH)- So, if these 
galaxies are truly low redshift galaxies, then their inferred star 
formation rates would be even lower. 

^ As we discuss in § 14.41 some of the submillimeter flux may be 
coming from accretion activity 
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Table 4. Properties of Submillimeter Galaxies with /s's-band Data 
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''Radio coordinates are f rom T able 2 of lChaoman et al.li2()05h . Five of these sources are also in lWang et"an 120041^ . The last source listed 
is only from Wang et ^ 120041) . 

''S'ssOMm fluxes are from Table 2 ofipiapman et al.' f2005'). Two of the sources are measured by^^^^^^^^^^^^^^ ^^^^ to have S'gso/im '"^ 
4.6 — 4.7 mjy, and are measured bv .Borvs et al. (2003) and Wa ng ct al. (2004) to have Sgso^im > 5 mjy, and for fairness we include these 
in the table. 

•^Spectroscopic redshift from lChapman et alj 1200511 . 

"^Soft and hard band fluxes are from 'Alexander et alj 120031) . 

"^This field indicates whether the submillimeter source satisfies the BX/BM, BzK /SF, and DRG selection criteria. 

^Inferred bolometric luminosity from iChaoman et al.l 120051'! . The mean bolometric luminosity of the 6 submillimeter galaxies with 
spectroscopic re dshifts 1.4 < z < 2.6 is (Lboi) ~ 9 X 10^^ Lq. 
SRedshift from IS winbank et al.l <20()1) . 



criteria designed to select star-forming galaxies. 

4.2. Passively Evolving Galaxies at z ^ 2 
Near-IR Colors 

We now turn to galaxies in our samples that ap- 
pear to have little or no current star formation. DRGs 
have SFRs that are comparable to those of BzK /SF 
and BX/BM galaxies with similar near-IR colors for 
(z — A')ab < 3. However, stacking the 13 DRGs with 
(z — K)ab > 3 results in a non-detection with an up- 
per limit of ~ 50 Mq yr~^ (Figure Stacking the 
X-ray emission from the 17 BzK /PE galaxies shows 
a similar turnover in the inferred average SFR around 
{z-K)ab ^ 3 (FigureUSl). The BzK colors of the DRGs 
(most of which only have i?-band limits) lie in the BzK 
color space expected for passively-evolving galaxies (Fig- 
ure I12|l . The stacking analysis confirms that these red 
DRGs and BzK/PFi galaxies have little current star for- 
mation compared with DRGs and _BzA'/PE galaxies with 
bluer near-IR colors (Figure I16|l. A similar X-ray stack- 
ing analysis bv lBrusa et al.1 1)2002|) yields no detection for 
passive EROs in the K20 survey from the Chandra Deep 
Field South data. 

The average J — K^ color of the 13 passively evolving 
DRGs (with iz-K)AB > 3) is {J-K^} = 2.98±0.59, and 
is comparable to the average J — Kg color of DRGs with 
bluer (z— if )ab colors, implying that the J-band is either 
in close proximity to the Balmer and 4000 A breaks or 
the band encompasses the breaks completely. This will 
occur for galaxies with a mean redshift (z) ~ 1.88 — 2.38 
(Figure OJ- In these cases, the (z — K)ab color will be 
more effective than the J— Kg color in culling those galax- 
ies with developed spectral breaks. The fraction of pas- 
sively evolving DRGs inferred from their lack of X-ray 
emission is 13/54 ^ 24%, which is in reasonable agree- 
ment to the passively evolving DRG fraction of ~ 30% 



found bv lForster Schreiber et all 1)20041) and lLabbe et alJ 
ll2005|) . 

Alternatively, galaxies with 2 < (z — K)ab < 3 must 
still be forming stars at a prodigious rate, as indicated 
by their stacked X-ray flux. The similarity in average 
Kg magnitude between galaxies with 2 < (z — K)ab < 3 
and {z — K) AB > 3 fFigure[TB|) suggests they have similar 
masses, and the difference in (z — K)ab color between 
the two samples simply reflects the presence of some rel- 
atively unobscured star formation in those galaxies with 
2 < (z - K)ab < 3. 

Stellar Populations 

FigurelT^furthcr illustrates the differences between the 
star-forming and passively evolving galaxies in terms of 
some physical models. The left panel of Figure [T51 shows 
the (z — K)ab versus {K — TO3.6/4.5/im)AB colors (near- 
IR/IRAC color diagram) for all galaxies, excluding di- 
rect X-ray detections ®. Because the SEDs of (non-AGN) 
galaxies considered here are expected to be relatively flat 
in fi, across the IRAC bands, we used 3.6 /im IRAC AB 
magnitudes for all sources that were not covered by the 
4.5 /im imaging. Al so shown in Fi gure IT^ are synthetic 
colors for [Bruzual fc Chariot ()2003) spectral templates 
at the mean redshifts of the BM ((z) ^ 1.7) and BX 
((z) ~ 2.2) samples, assuming constant star formation, 
E{B -V) = and 0.3, and a fGaTzetti erall ^00) red- 
dening law. The bulk of BX/BM and BzK/SF galax- 
ies generally fall within the region of color space ex- 
pected for constant star forming galaxies with moderate 
extinction of E{B — V) ^ 0.15 and ages of ~ 1 Gyr. 
These values are consistent with those derived from de- 
tailed spectral modeling of BX/BM galaxies and LBGs 

^ ILabbe eTall <2(X)5() propose a similar diagram to differentiate 
DRGs from other (e.g., star- forming) populations of z ~ 2 galaxies. 
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bv iShaplev et al.1 1)2005(1 . Much of the scatter of star- 
forming galaxies to the left and right of the CSF models 
for E{B — V) = 0.15 — 0.30 is a result of photometric un- 
certainty, particularly in the {K — TO3.6/4.5/^m)AB color, 
since we include galaxies with formal IRAC uncertainties 
of 0.5 mag. In addition, some of the scatter for objects 
with blue {z — K)ab colors arises from interlopers. The 
more interesting aspect of FigurellSb is that the constant 
star formation models cannot account for the colors of 
objects with (z — K)ab ^ 3: these objects must have 
ages less than the age of the universe at z ~ 2 (~ 3 Gyr) 
and simultaneously have modest E{B — V) — and hence 
modest current SFRs< 190 M© yr~-^ — such that they 
remain undetected as soft-band X-ray sources. The im- 
portant result is that, similar to SED fitting, the X-ray 
stacking analysis allows us to rule out certain star forma- 
tion histories. For the PE galaxies considered here, the 
X-ray data rule out the constant star formation models. 
The benefit of X-ray data is that we can quantify the cur- 
rent SFR independent of extinction and the degeneracies 
which plague SED fitting. 

The only models that can account for the colors of ob- 
jects with (z — K)ab ^ 3 are those with declining star 
formation histories (right panel of Figure I18|) . For ex- 
ample, DRG/PEs at z ^ 2.2 have colors that can be 
reproduced by dust free models {E{B — V) = 0.0) with 
star formation decay timescales between r = 10 Myr 
(instantaneous burst) and r ~ 700 Myr^. While the up- 
per limit on the current average SFRs of DRG/PEs of 
50 Mq yr~^ (Figure 116)1 does not help us further con- 
strain the star formation history to a narrower range in 
r, the fact that DRG/PEs are still detected at z-band 
suggests that single stellar population models with small 
r are unrealistic (e.g., ages greater than 1 Gyr imply 
> 100 e-folding times for the instantaneous burst model 
making such an object undetected at z-band). If there 
is ongoing low level star formation activity, then a two 
component model with a underlying old stellar popula- 
tion and a recent star formation episode may be required 
Ce-g.- IYan eHIMoO^ . 

The red (z — K)ab > 3 colors of BzK/PFi galaxies 
can be reproduced by models with r ^ 10 — 300 Myr 
(for T much larger than 300 Myr, the models over pre- 
dict the current star formation rate). The X-ray data 
indicate that these BzK /PE galaxies have an average 
current SFR of ~ 28 M© yr"!. For an age of 1 Gyr 
and T — 300 Myr, this implies an "initial" SFR of 
^ 800 M© yr^i at z 2.4. This initial SFR is compara- 
ble to that of SMGs and the implied formation redshift 
is close to the median redshift of SMGs (z ~ 2.2), sug- 
gesting SMGs could be plausible progenitors of BzK/PE 
galaxies if the single component model correctly de- 
scribed the star form ation histories of these galaxies (e.g., 
ICimatti et al.l I2004D . The formation redshifts can be 
pushed back in time to significantly earlier epochs z > 3.5 
if one assumes a more physically motivated "truncated" 
star formation history that models the effects of feedback 
in halting star formation (Daddi et al. 2005,) . Nonethe- 
less, the simplistic example above illustrates how X-ray 
estimates of the bolometric SFRs of galaxies can be com- 
bined with the results of stellar population modeling to 

® We rule out ages that are greater than the age of the universe 
at 2 ~ 2.2, giving an upper limit on the age of ~ 3 Gyr. 



indicate the likely progenitors of such galaxies. In sum- 
mary, we have demonstrated the utility of stacked X-ray 
data as a powerful constraint on the results from stel- 
lar population modeling. The X-ray data indicate that 
galaxies with (z — K)ab ^ 3 have SEDs that are consis- 
tent with declining star formation history models. Other 
studies also show that the SEDs of BzK/PE galaxies 
can be adequately described by declining star formation 
history models (e.g., Daddi ct al. 2005). The deep X-ray 
data confirm these results and further allow us to con- 
strain the current SFRs of passively evolving and star- 
forming galaxies independent of the degeneracies associ- 
ated with stellar population modeling. 

Space Densities 

To conclude this section, we note that objects se- 
lected by their J — Kg colors appear to include a sub- 
stantial number of passively evolving galaxies at red- 
shifts z > 2 (although, as we pointed out above, the 
(z — K)ab color may be a more effective means of de- 
termining which DRGs are passively evolving). The 
space density implied by the 4 BzK/PFi galaxies with 
(z - K)ab > 3 and ifs < 21 is ~ 3 X 10"^ Mpc'^, as- 
suming a boxcar (or top-hat) selection function between 
redshifts 1.4 < z < 2.0 and an area of ~ 72.3 arcmin^. If 
we include all 17 BzK/PE galaxies (i.e., including those 
with {z — K)ab < 3 and those four that are directly- 
detected in X-rays) with Kg < 21, we find a space den- 
sity of 1.3 X 10^"* Mpc^'^. Given the st rong cluster- 
ing o bserving for BzE/PE galaxies (e.g., iDaddi et alJ 
12005)) . our estimate is in good agreement with the value 
of 1.8 X lO"'* Mpc-3 obtained for the Daddi ct al. { 206^ 
sample of bona- fide BzK/PE objects in the Hubble Ultra 
Deep Field (UDF), an area which is 6 times smaller than 
the area considered in our analysis (i.e., 12.2 arcmin^ in 
^addi et al. ( 2005) versus 72.3 arcmin considered here). 

All of the i?zi^/PE galaxies would have been detected 
with Ks < 21 and {z—K)ab > 3, assuming the PE model 
shown in Figure 118b . if they were at the mean redshift 
assumed for DRG/PEs of (z) — 2.2. The space density 
of the 13 DRG/PEs with (z - K)ab > 3 is - 9 x 10"^ 
to Ks = 21, assuming a boxcar selection function be- 
tween between redshifts 2.0 < z < 2.6 and an area of 

72.3 arcmin^. Here, we have assumed that the range 
in redshifts of DRG/PEs (2.0 < z < 2.6) is similar to 
that of all DRG galaxies based on the spectroscopic red- 
shift distribution shown in Figure ^ Our estimate of the 
D RG/PE space densit y is comparable to that obtained 
bv iLabbe et al.l l)2005|) after restricting their sample to 
Kg < 21 (yielding 1 object over ~ 5 arcmin^), and as- 
suming a volume between redshifts 2.0 < z < 2.6). 

It is worth noting that the BzK /PE and DRG/PE 
populations appear to be highly clustered (e.g., 
IDaddi et alJl2005llvan Dokkum et aTll2004)) and this will 
likely affect their density estimates over small volumes. 
While the space densities derived here are in rough agree- 
ment with other studies, our estimates have been derived 
over a much larger volume (by a factor of 6 — 14) than 
any previous study and will be less susceptible to varia- 
tions in density due to clustering. We caution the reader 
that the density estimates may still be uncertain given 
that the redshift distribution of DRGs with z > 2.6 is 
not well sampled, even in the large spectroscopic dataset 
of DRGs considered here. 
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Fig. 18. — (a) (z - K)ab versus (i<' - rris g/4 5^jn)AB colors for Ks < 21 galaxies in the BX/BM, Bzft'/SF, BzK/PE, and DRG samples. 
Also shown are lBruzual Hr, C hariot" ( 200?) spectral templates for ages 1 Myr to 3 Gyr at the mean rcdshifts of the BX {(z) 2.2) and BM 
({2) ~ 1.7) samples, assuming constant star formation, E{B — V) = 0.0 and 0.3, and the|Calzetti ct al. 12000) reddening law. The models 
assume a Salpeter IMF from 0.1 — 100 Mq and solar metallicity. The photometric scatter is large given that we include galaxies with IRAC 
uncertainties up to 0.5 mag, and this accounts for the large spread in the {K — "13 g/4 5^jn)AB colors of the star-forming candidates; (b) 
same as (a) except we show the spectral templates for a model with r = 200 Myr and E{B — V) = 0.0. The shaded region selects IRAC 
Extremely Red Objects (lERGsi .Yan ct al...2004) . 



Taken at face value, the density estimates derived 
above suggest a significant presence of passively evolving 
Ks < 21 gal axies at redshifts z > 2. This result contrasts 
with that of lDaddi et alJ l)2005|) who argue for a rapid de- 
crease in the number density of passively evolving galax- 
ies at redshifts z > 2 based on BzK/PFi selection. The 
dropoff in space density of passively ev olving galaxie s 
at redshifts z > 2, as sugg ested by Daddi et alJ 1|2005l) . 
may be an artifact of the BzK/PFi selection function 
which, based on previously published redshift distribu- 
tions and shown in Figure Q appears to miss passively 
evolving galaxies at redshifts z > 2, even when using the 
very deep i3-band of this study. Figure IT^ shows that all 
of the DRGs that cluster around the BzK/PFi selection 
window have limits in i?-band and very few actually have 
limits that would for certain place them in the BzK/PE 
window. Photometric scatter for those galaxies with very 
faint (or no) -B-band detections is likely to be significant 
for these galaxies. These results suggest that the depth 
of the B-band data is the determining factor in whether 
BzK/PE selection culls galaxies with redshifts 2; > 2 or 
not. Because the depth of the photometry is an issue for 
the BzK/PE selection, it becomes difficult to accurately 
quantifying with a single selection criteria the dropoff 
in space density of passively evolving galaxies between 
z < 2 and z > 2. We can avoid the need for excessively 
deep B-band data to select passively evolving galaxies 
with redshifts z > 2 by simply selecting them using a sin- 
gle color, J— Kg or {z~K)ab- The stacked X-ray results 
show that a subsample of DRGs have very little star for- 
mation, suggesting that passively evolving galaxies have 
a significant presence at epochs earlier than z = 2. The 
inferred ages of DRGs would imply formation redshifts 
of z - 5 (Lab be et al.l l2005). 

We conclude this section by noting that there are sev- 



eral ga laxies in the HS1700-t-643 sample of lShaplev et alJ 
((20Q3), and r nany more in op tically selected samples in 
general Ce.g.. lErb et alJl200al . which have old ages and 
early formation redshifts similar to those of the passively 
evolving BzK/PE and DRG/PE galaxies discussed here. 
In order to reproduce the observed SEDs for such ob- 
jects, the current SFR must be much smaller (but still 
detectable in the case of optically selected galaxies) than 
the past average SFR. 

4.3. Selecting Massive Galaxies 

As discussed above, DRGs with (z — K)ab > 3 ap- 
pear to be passively evolving based on their (lack of) 
stacked X-ray flux and their colors with respect to mod- 
els with declining star formation histories. The X-ray 
data indicate that BzK/PE galaxies also appear to be 
well described by declining star formation h istories, con- 
sistent with the SED modeling results of iDaddi et alJ 
l)2005ft . The stellar mass estimates of these PE galax- 
ies will be presented elsewhere. Here, we simply mention 
that several existing studies of the stellar populations 
of BzK/PE and DRG galaxies with Kj < 20 indicate 
they have masses > 10" (e.g., iDaddi et al. 200i 
iForster Schreiber et al.l l2004|) . In addition, Yan et aQ 
()2004|) recentlv analyzed the stellar populations of IRAC- 
selected Extremely Red Objects (lEROs), selected to 
have /i/(3.6/im)//y(z85o) > 20 (or, equivalently, (z — 
3.6/im)AB > 3.25). Spectral modeling indicates these 
sources lie at redshifts 1.6 < z < 2.9, are relatively old 
(1.5 — 3.5 Gyr), and require an evolved stellar popula- 
tion to fit the observed SEDs. Almost all of the PE 
galaxies with (z — K)ab > 3 satisfy the lERO criteria 
(shaded region of FigurclTM)). Furthermore, the 7?,-band 
detections and limits for PEs with (z — /^)ab > 3 imply 
TZ ~ Ks> 5.3, satisfying the ERO criteria. 
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The inferred stellar masses of the lShaplev eTall lEM) 
sample of optically-selected galaxies in the HS1700-I-643 
field are shown in Figure El For compariso n, we also 
show the inferred stellar masses from the lYan et al. 
(1200411 sa mple of lEROs^". Both the iShaplev et al 
|2OO90 and lYan ~al. (2004) samples take advantage of 
the longer wavelength IRAC data to constrain the stel- 
lar masses, and the typical uncertainty in mass is ^ 40% 
for objects with M* < 10" Mq and < 20% for objects 
with M* > 10" M0. The lERO stellar masses have 
been multiplied by 1.7 to convert from a Chabrier to 
Salpeter IMF. The scatter in Figure 111?! reflects the large 
(> 1 magnitude) variation in the mass-to-light (M/L) ra- 
tio for objects with a given rest-frame optical luminosity 
in the BX and ZERO samples. For magnitudes brighter 
than our DRG completeness limit of Ks = 21, ~ 16% 
of BX galaxies have stellar masses M* > 10^^ Mq. For 
Ks < 20 BX galaxies, the fraction with M* > 10" Mq 
is ~ 40%. So, while optically-selected BX galaxies have 
a lower mean stellar mass than lEROs, there is certainly 
a subsample of BX galaxies which have masses compara- 
ble to the most massive lERO galaxies. We note that the 
stellar mass distributions of > 21 BXs and lEROs do 
not overlap: the TZ — Kg colors of lEROs are too red for 
them to be included in the optically-selected sample. 
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Fig. 19. — Infe rred stellar masses of BX and "MD" objects 
from the lijhaplev e t al. 12005) sample with spec troscopic redshifts 
1.5 < Zspec < 2.9 (circles), and lEROs from the lYan et'ai] l200l) 
sample with photometric redshifts 1.6 < Zphot < 2.9 (squares, and 
triangles for those with Ks-hand limits). The stellar masses from 
lYan et al. 12004) have been multiplied by 1.7 to convert from a 
Chabrier to Salpeter IMF. The dashed vertical line denotes the 
limit brighter than which we are complete for lEROs (i.e., those 
DRGs with (z — -/<')ab > 3). A subset of BX galaxies have stellar 
masses similar to those of the lERO sample. The scatter in stellar 
masses reflects at least a magnitude variation in the M/L ratio of 
BX/MD and lERO objects at a given rest- frame optical luminosity. 

The X-ray stacking results indicate that BX/BM 
galaxies with Kg < 20 have prodigious star formation 
rates, while lEROs (i.e., those DRGs with (z—K) ab > 3) 

H. Yan, private communication 



have very little star formation. Therefore, a simple inter- 
pretation is that optical surveys include objects that are 
as massive (M* > 10^^ Mq) as those selected in near- 
IR surveys, with the only requirement that the galax- 
ies have some unobscured star formation. The range of 
star formation rates (uncorrected for extinction) found 
for BX/BM galaxies is 3 — 60 Mq yr~^, and it is likely 
that massive galaxies with at least a little unobscured 
star formation can be BX/BM selected. This may be 
the only significant difference between optical and near- 
IR selected massive galaxies, and the difference in SFR 
may be temporal. These criteria typically fail to select 
passively evolving galaxies at z ~ 2 as they have al- 
ready settled to a quiescent stage. This does not mean 
that such massive galaxies will never appear in optical 
surveys. For example, a subsequent accretion event at 
z < 2 could elevate the star formation activity in an oth- 
erwise passively evolving massive galaxy, thus bringing 
it into the optical sample. Nonetheless, the DRG and 
BzK /PE criteria add to the census of galaxies at z ~ 2 
by selecting passively evolving galaxies that have stellar 
masses similar to the most massive galaxies selected in 
the rest-frame UV. 

4.4. Star Formation Rate Density at z ^ 2 

Contribution from Optical and Near-IR Selected 
Samples 

We can roughly estimate the contribution of BX/BM, 
BzK/SF, and DRG/SF galaxies to the extinction-free 
star formation rate density (SFRD) at z ~ 2^^. The 
BzK /SF criteria are designed to select galaxies with 
redshifts 1.4 < z < 2.6. The similarity in surface 
densities, volumes probed, and SFRs of galaxies in the 
BzK /SF and BX/BM samples implies that their contri- 
bution to the SFRD will be comparable for objects with 
Ks < 22. The redshift distribution of A'^ < 21 DRGs 
from within our own sample is reasonably well defined 
over this redshift range (cf., Figure so we can esti- 
mate the added contribution of 80% of the DRGs with 
redshifts 2.0 < z < 2.6 to the SFRD between redshifts 
1.4 < z < 2.6. Figure 1^ and Tabled show the cumu- 
lative contribution to the SFRD of BzK /SF, BX/BM, 
and DRG/SF galaxies. The points in Figure QUI are not 
independent of each other due to the overlap between the 
samples (e.g.. Figure^}. Also shown in Figurel^bv the 
shaded region is the inferred total SFRD assuming the 
overlap fractions of Figure ^1 and counting all objects 
once. The results indicate that BX/BM selection would 
miss approximately one-third of the total SFRD from 
BX/BM and BzK /SF galaxies to K^ = 22 and DRG/SF 
Kg < 21 galaxies combined. We remind the reader that 
much of the incompleteness of the BX/BM sample with 
respect to that of the BzK /SF sample (and vice versa) 
results from photometric scattering (e.g., Figures lTTl and 
I12|) . Monte Carlo simulations can be used to quantify 
the biases of such photometri c inaccuracy and th us cor- 
rect for incompleteness (e.g., IReddv et al.ll20()5|) . The 
total SFRD in the interval 1.4 < z < 2.6 for BX/BM 

Although the BzK /PE galaxies do have detectable X-ray 
emission (e.g., Figure lTSt . their contribution to the SFRD is mini- 
mal given that their space density is a factor of 5 smaller than that 
of the BzK /SF and BX/BM galaxies to Ks = 21 in the redshift 
range 1.4 < z < 2.0. 
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Table 5. Cumulative Contributions to the SFRD Between 1A < z < 2.6 





Range 




Sample^ 




1-1 


P 

(arcmin } 


QT7''D T~ld 

(Mq yr"^ Mpc~^) 


K 


< 20.0 


RX /RM 






1.4-2.6 


0.44 


\J .\J \-\J _1_ W.iJL/O 












1.4-2.6 


0.62 


n ni a + n nr)4 






DRG/SF 






2.0-2.6 


0.15 


0.008 lb 0.002 






BzK/S¥- 


-BX/BM 




1.4-2.6 


0.35 


0.007 lb 0.002 










RY /RM 


2.0—2.6 


08 


n 004 -1- 001 






Total 






2 0—2 6 


U.o / 


n 097 -1- OOfl 


K= 


< 20.5 


RX /RM 






1.4-2.6 


1.13 


n fiQA _|_ n nnn 

VJ.L/Ort _1_ L/.VJL/J 












1.4-2.6 


1.35 


n 032 + n06 






DRG/SF 






2.0-2.6 


0.30 


0.013 lb 0.003 






BzKj'SF- 


-BX/BM 




1.4-2.6 


0.55 


0.012 lb 0.002 






DRG /SF— 


-R7 /SF- 


RY /RM 


2.0—2.6 


0.25 


n OOQ -1- 009 






Total 






2 0—2 6 


1 0*^ 
1 . yo 


n o^^f^ -1- 01 1 




< 21.0 


BX/BM 






1.4-2.6 


2.15 


0.045 lb 0.010 






B^A'/SF 






1.4-2.6 


2.34 


0.044 lb 0.009 






DRn /SF 






2.0—2.6 


0.58 


n 090 -1- 004 






DZl\ / Ol^ — 






1 4—2 6 


n 74 


01 A -1- 009 






DRC/SF— 


J—} J.\. 1 kJ J- 


— RX /RM 


2.0-2.6 


0.41 


01 H- 009 

VJ.W-LvJ _1_ \J .\J\J £j 






Total 






2.0-2.6 


3.30 


0.072 ±0.014 


Ks 


< 21.5 


BX/BM 






1.4-2.6 


3.17 


0.057 ±0.013 






Sz/sT/SF 






1.4-2.6 


3.61 


0.052 ±0.010 






Bz/sT/SF- 


-BX/BM 




1.4-2.6 


1.71 


0.017 ±0.002 






TotaP 






2.0-2.6 


5.29 


0.087 ±0.017 




< 22.0 


BX/BM 






1.4-2.6 


4.31 


0.069 ±0.015 






BzX/SF 






1.4-2.6 


4.93 


0.068 ±0.014 






BzA'/SF- 


-BX/BM 




1.4-2.6 


2.55 


0.020 ± 0.003 






TotaF 






2.0-2.6 


7.27 


0.102 ±0.021 



^The Bzii'/SF — BX/BM sample represents the set of objects that are Szi^/SF-selected, but not BX/BM-selected. Similarly, the 
DRG/SF — BzKI'SP — BX/BM sample represents the set of DRGs that are not selected by either the BziC/SF or BX/BM criteria. 
'^Redshift range of sample. 

•^Surface density of photometric objects after removing interlopers and directly detected X-ray sources that are likely AGN. The number 
of BX/BM objects is calculated assuming the spectroscopic and interloper fractions from Table The overlap fractions are taken from 
Figure [Tni We assume a field area of ~ 72.3 arcmin^ to compute surface densities. 
Assumes the average SFRs shown in Figure [T^ 

<=This includes the contribution from the DRG/SF— BzA'/SF— BX/BM sample for < 21. 



and BzK /SF galaxies to = 22 and DRG/SF galax- 
ies to Kg = 21, taking into account the overlap between 
the samples, is ~ 0.10 ± 0.02 Mq yr^^ Mpc^^. Approx- 
imately 30% of this comes from galaxies with Kg < 20 
(Table [SJ. Optically-selected galaxies to 7^ = 25.5 and 
Ks = 22.0 and BzK/SF galaxies to = 22.0 (with 
significant overlap between the two samples) account 
for - 87% of the total SFRD quoted above. DRGs to 
Kg = 21 that are not selected by the BX/BM or BzK /SF 
criteria contribute the remaining ^ 13%. We note that 
the number ~ 0.10 ± 0.02 Mq yr'^ Mpc"^ does not 
include the 5 radio-selected SMGs to Sssoiim ^ 5 mJy 
that are near-IR and/or optically-selected since we re- 
moved the directly-detected hard-band X-ray sources in 
computing the SFRD. If we add these 5 radio-selected 
SMGs that are present in the optical and near-IR sam- 
ples (all of which have Kg < 21), then the total SFRD 
contributed by the BX/BM and BzK/SF objects to 
Kg = 22.0 and DRG/SF galaxies to Kg = 21.0 is 
0.15 ± 0.03 Mq yr^i Mpc~^ (see next section). 

Contribution from Radio- Selected Subniillimeter 
Galaxies 

We conclude this section by briefly considering the 
contribution of radio-selected submillimeter galaxies 



(SMGs) with S'ssoA.m > 5 mJy to the SFRD. All but 
one of the radio-selected SMGs summarized in Table 0] 
are directly detected in either the soft or hard band 
and are likely associated with AGN. Stacking the X-ray 
emission for the 5 radio-selected SMGs with redshifts 
1.4 < z < 2.6 in Table H yields an average inferred 
SFR of ^ 2900 Mq yr~\ and this value should be re- 
garded as an upper limit given that the X-ray emission 
is likely contaminated by AGN. On the other hand, the 
average bolometric luminosity of the 5 SMGs, as derived 
from their submillimeter flux, is (iboi) ~ 9 x 10^^ Lq. 
If we assume th at 30 % of Lhoi arises from AGN (e.g., 
'Ghap man et alJ |2005l. lAlexander et a l. 2005"l. then the 
implied SFR is ~ 1000 Mq yr"^ If we take at face 
value the assertion that 30% of the bolometric luminos- 
ity of submillimeter galaxies comes from AGN, then this 
means the X-ray emission would overestimate the aver- 
age SFR of SMGs by a factor of ^ 3. In other words, 
only one-third of the X-ray emission from SMGs would 
result from star-formation, and the remaining two-thirds 
would come from AGN. 

To determine the additional SFRD provided by radio- 
selected SMGs with 5850pm ^ 5 mJy, we must account 
for their overlap with the optical and near-IR samples. 
The data in TableHshow that there are 4 of 9 SMGs that 
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Fig. 20. — Cumulative star formation rate density (SFRD) as a 
function of Kb magnitude for BX/BM (circles), BzK/SF (trian- 
gles), and DRG/SF (squares) galaxies with redshifts 1.4 < z < 2.6. 
The points are not independent of each other given the overlap 
between the samples. The shaded region denotes the total cumu- 
lative SFRD when counting overlap objects once. The total SFRD 
to Ks = 22 includes DRGs with Kg < 21. The error bars reflect 
the Poisson error and uncertainty in star formation rate added in 
quadrature, but do not reflect systematic errors associated with, 
e.g., photometric scattering. 



are not selected by the optical and/or near-IR criteria. 
All 4 of these galaxies have relatively low redshifts z < 1 
and will obviously not contribute to the SFRD between 
redshifts 1.4 < z < 2.6. Alternatively, of the 5 SMGs 
that are spectroscopically confirmed to lie at redshifts 
1.4 < z < 2.6, all are selected by either the BX/BM, 
BzK/SF, or DRG criteria (and sometimes by more than 
one set of criteria). 

Because of the non-uniform coverage of the submil- 
limeter observations, we must rely on the published sub- 
millimeter number counts to estimate the effective sur- 
face density probed by the 9 radio-selected SMGs with 
•S'sso/iin ^ 5 mJy hsted in Table IH Accordi ng to the mod- 
els shown in Figure 4 of lSmailet af]ll200l . we should ex- 



pect to find ~ 0.25 sources arcmin" to S'sso/^m ~ 5 mJy. 
Neglecting cosmic variance, the nine observed radio- 
selected SMGs to 5850^111 ~ 5 mJy then imply an effec- 
tive surface area of ~ 36 arc min^ . The spectroscopi c red- 
shifts compiled in Table 2 of lChanman et all l)2005D indi- 
cate that « 44% of the radio-selected SMGs to S'ssoAtm ~ 
5 mJy lie at redshifts outside the range 1.4 < z < 2.6. 
If we assume a Poisson distribution of sources, then the 
total number of SMGs to S'sso^m ^ 5 mJy could be as 
high 9 + \/9 =12. If we assume that the fraction of 
interlopers among the 3 unobserved objects is similar 
to the fraction of interlopers among objects which are 
observed, then we expect an upper limit of two SMGs 
with S'sso/iin ~ 5 mJy that are unobserved and that lie 
between redshifts 1.4 < z < 2.6. If we conservatively 
assume that these two sources are not selected by the 
optical and/or near-IR criteria, and they have bolomet- 
ric SFRs of 1500 Mq yr~i (similar to the average SFR 



found for the 5 spectroscopically confirmed radio-selected 
SMGs in Table then the inferred additional SFRD 
provided by these two SMGs would be ~ 3000 Mq yr~^ 
divided by the volume subtended by 36 arcmin^ at red- 
shifts lA < z < 2.6, or ~ 0.022 Mq yr"! Mpc~3_ 
We note that this should be treated as an upper limit 
for several reasons. First, we have assumed the maxi- 
mum number of unobserved sources allowed by Poisson 
statistics. Second, we have assumed an interloper frac- 
tion among these unobserved sources that is the same 
for the observed sources. In general, one might expect 
the contamination fraction to be higher among the gen- 
eral SMG population to S'sso/im ~ 5 mJy (where an ac- 
curate radio position may not be known) than would 
be inferred from the radio-selected SMG surveys. Fi- 
nally, we have assumed that all of the unobserved sources 
cannot be selected by their optical and/or near-IR col- 
ors. Neglecting any overlap, radio-selected SMGs to 
•S'ssoAim ^ 5 mJy contribute ~ 0.05 Mq yr~^ Mpc^^ to 
the SFRD between redshifts 1.4 < z < 2.6. However, 
our conservative calculation indicates that radio-selected 
SMGs to SsbOfim ^ 5 mJy that are not selected by opti- 
cal (BX/BM) and/or near-IR {BzK and/or DRG) sur- 
veys make a small (< 0.022 Mq yr-^ Mpc"^ or < 15%) 
additional contribution to the SFRD between redshifts 
1.4 < z < 2.6. 

5. CONCLUSIONS 

We have taken advantage of the extensive multi- 
wavelength data in the GOODS-North field to select 
galaxies at z ^ 2 based on their optical and near-IR 
colors and to compare them in a consistent manner. Our 
own ground-based optical and near-IR images are used 
to select galaxies based on their UnGTZ, BzK, and J—Ks 
colors. Additional rest-frame UV spectroscopy for 25% 
of optically selected candidates allows us to quantify the 
redshift selection functions for the various samples. We 
use the deep Chandra 2 Ms X-ray data to determine 
the influence of AGN and estimate bolometric SFRs for 
galaxies in the optical and near-IR samples. We also 
use the deep Spitzer/IRAC data in the GOODS-North 
field in considering the stellar populations and masses of 
galaxies selected in different samples. Our analysis em- 
ploys the same multi-wavelength data for a single field 
(GOODS-North), using the same photometric measure- 
ment techniques, for a consistent comparison between 
galaxies selected by their optical and near-IR colors. Our 
main conclusions are as follows: 

1. Star- forming galaxies at z ^ 2 selected by their 
UnGTZ colors (i.e., BX/BM galaxies) and their BzK 
colors (i.e., BzK/SF galaxies) have optical and near- 
IR color distributions that indicate significant overlap 
(~ 70 — 80%) between the two samples. Photomet- 
ric scatter could account for the colors of at least half 
of those galaxies missed by one set of criteria or the 
other. The BzK/SF criteria are less efficient in selecting 
(younger) Kg > 21 galaxies at redshifts 1.4 < z < 2.6, 
while the BX/BM criteria are less efficient in selecting 
near-IR bright (e.g., Kg < 20) objects. Distant red galax- 
ies (DRGs; including both reddened star-forming and 
passively evolving galaxies) selected to have J — Kg > 2.3 
show near-IR colors that are 1 — 1.5 magnitudes redder 
than for samples of star-forming galaxies. Criteria aimed 
at selecting passively evolving galaxies based on their 
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BzK colors (i.e., BzK /PE galaxies) by design have red 
near-IR colors, but we find that the redshift distributions 
of BzK/PE galaxies and DRGs have very little overlap. 

2. The deep X-ray data show that almost all of the 
directly detected X-ray sources in the samples have hard 
band emission and X-ray/optical flux ratios indicating 
they are likely AGN. Much of this AGN contamination 
occurs for magnitudes Kg < 20. We identify 5 objects 
that are detected in the soft band X-ray data and are 
likely star-forming galaxies based on their absence in the 
hard-band X-ray data, optical magnitudes i? > 22, and 
absence of obvious AGN features for those with spec- 
tra. We stacked the X-ray data for all likely star-forming 
galaxies (i.e., those undetected in X-rays and the 5 galax- 
ies discussed above), excluding likely AGN. The stacking 
analysis shows that the star formation rate (SFR) dis- 
tributions of BX/BM and BzK /SF galaxies and DRGs 
are very similar as a function of Kg magnitude. Galaxies 
with < 20 have average SFRs of - 120 Mq yr-\ a 
factor of two to three higher than Ks > 20.5 galaxies. 
Previous studies point to a similarity in the metallici- 
ties, clustering, and stellar masses of Ki, < 20 optical 
and near-IR selec ted galaxies (e.s.. iShaolev et ;il]l2fl04 
lAdelberger et alJ l20051. In this work we show that the 
bolometric SFRs of optical and near-IR selected galax- 
ies are also very similar when subjected to a common 
near-IR magnitude. 

3. Near-IR selection of star forming galaxies should be 
more immune to the effects of dust obscuration than opti- 
cal surveys. However, the BX/BM, BzK /SF, and DRG 
samples show very similar SFRs as a function of near- 
IR color for galaxies with (z — i^)AB < 3. The SFRs 
inferred for BzK /SF galaxies which are not optically- 
selected are very similar to BzK /SF galaxies which do 
satisfy the optical criteria, suggesting that star-forming 
galaxies in near-IR samples that are missed by opti- 
cal criteria do not harbor large numbers of heavily red- 
dened galaxies. Furthermore, the optical and BzK /SF 
samples host an approximately equal number of submil- 
limeter galaxies (SMGs). Despite their very high SFRs 
(~ 1000 M0 yr~^), SMGs missed by optical and near- 
IR selection have insufficient space densities to make a 
significant contribution (< 10%, see conclusion 6 be- 
low) to the total census of star formation at redshifts 
1.4 < z < 2.6. 

4. We identify a population of extremely red BzK 
and DRG galaxies with (z — K)ab ^ 3. The stacked 
X-ray data indicate these red galaxies have little, if any, 
current star formation. The absence of X-ray emission 
from these objects also suggests that low luminosity AGN 
and low mass X-ray binaries contribute little X-ray emis- 
sion in star-forming galaxies compared with the emis- 
sion produced from more direct tracers of the current 
star formation rate, such as high mass X-ray binaries. 



We further demonstrate the utility of deep X-ray data 
to constrain the stellar populations of these extremely 
red galaxies, and find that they must be described by 
declining star formation histories. Almost all of these 
passively e volving galaxies satisfy the lERO criteria of 
lYan et all (f2004.) . We find that optical selection includes 
a subset of galaxies with stellar masses similar to those 
inferred for lEROs, but which are forming stars at a 
prodigious rate . The stellar mass estimates from SED 
mode ling (e.g., lYan et a T 2004, For ster Schreiber et alJ 
l2004j) and bolometric SFR estimates from the deep X- 
ray data (this work) suggest that the presence or absence 
of star formation may be the only significant difference 
between optical and near-IR selected massive galaxies 
(M* > 10" Mq), and the difference in SFR may be 
temporal. 

5. We find evidence for a significant presence of pas- 
sively evolving galaxies at redshifts z > 2 compared with 
their space density at lower redshifts, 1.4 < z < 2.0. Our 
analysis suggests that a single color technique using the 
(z — K)ab or J — Kg color allows for a more practical 
method selecting passively evolving galaxies with z > 2 
than the BzK /PE criteria, as the latter would require ex- 
cessively deep i?-band data to accurately determine the 
space densities of passively evolving galaxies at z > 2. 

6. Finally, we consider the contribution of optical 
and near-IR selected galaxies to the SFRD at z ~ 2, 
taking into account the overlap between the samples 
and their respective redshift distributions. We find 
that BX/BM and BzK /SF galaxies to K^ = 22, and 
DRG galaxies to Kg = 21, account for an SFRD of 
- 0.10 ± 0.02 Mq yr-i Mpc'^ between redshifts 1.4 < 
z < 2.6. Approximately 87% of this total comes from 
optically-selected galaxies to TZ — 25.5 and Kg = 22 
and near-IR selected BzK galaxies to Kg = 22, and 
13% from Ks < 21 DRGs not selected by the BX/BM 
or BzK criteria. Of the known radio-selected SMGs to 
Ss50fj,m ^ 4 mJy in the GOODS-N field with redshifts 
1.4 < z < 2.6, > 80% could be selected by the BX/BM, 
BzK, and/or DRG criteria. 
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